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Abstract 

A  comprehensive  treatment  of  the  theory  of  rf  energy 
absorption  in  inhomogeneous  plasmas  by  means  of  Alfv£n  waves 
is  presented.   A  straight  cylindrical  plasma  column  whose 
motion  is  governed  by  linearized  ideal  magnetohydrodynamics 
(MHD)  is  assumed  to  be  surrounded  by  a  vacuum  region  in 
which  an  ideal  coil  consisting  of  surface  currents  is  placed. 
The  coil  is  oscillated  at  frequencies  which  lie  in  the  MHD 
continuous  spectrum  of  local  Alfv£n  frequencies  of  the  inhomo- 
geneous  plasma  leading  to  the  occurrence  of  local  Alfv£n 
resonances.  A  rigorous  self-consistent  analysis  of  the 
resulting  resonant  energy  absorption  is  presented  along 
with  calculations  relevant  to  present  and  future  fusion 
devices.   Specifically,  the  experimentally  useful  coil 
impedance  is  calculated  for  a  plasma-vacuum-coil-wall  system 
with  realistic  nonuniform  plasma  profiles  corresponding  to 
present  day  Tokamak  and  Scyllac  devices.   Impressively  large 
coupling  and  absorption  is  indicated  for  such  devices 
independent  of  the  value  of  total  plasma  8.   Significant 
differences  occur  however  for  variations  in  the  spatial  plasma 
profiles.   Consequently  accurate  scaling  estimates  of  Alfve'n 
wave  absorption  are  not  possible  without  solving  the  complete 
nonuniform   problem. 


1 .   Introduction 

One  of  the  most  important  problems  facing  the  development  of 
Tokamak  fusion  plasma  experiments  is  the  effective  plasma  heating 
over  and  above  that  obtainable  from  Joule  heating.   This  paper 
addresses  itself  to  this  problem.   In  particular,  calculations 
are  presented  for  a  new  heating  mechanism  which  may  be  more 
attractive  than  current  conventional  techniques. 

Of  the  many  candidates  for  Tokamak  heating  which  are  currently 
in  vogue  are  neutral  beam  injection  and  rf  plasma  heating  at  ion 
cyclotron  (ICRH)  and  lower  hybrid  frequencies  (LHRH) .  Neutral 
beam  injection  has  proven  effective  in  small  machines  such  as 
the  ATC  and  ORMAK  where  the  beam  only  produces  small  perturba- 
tions on  the  background  plasma  equilibrium.   Fusion  plasma 
experiments,  on  the  other  hand,  due  to  the  very  large  beam  energy 
and  currents,  may  be  significantly  disturbed  by  the  beam  such 
that  the  equilibrium  would  be  altered.   Ion  cyclotron  resonant 
heating  has  also  been  shown  to  work  in  the  ST  plasma.  The  rf 
power  densities  in  the  ST  experiment  were  of   the  order  of  100  kw, 
and  significant  ion  heating  was  observed.   Extrapolating  again  to 
fusion  plasma  parameters,  it  is  questionable  whether  such  large  rf 
energies  can  be  generated  at  the  necessary  high  frequencies  .    From 
the  standpoint  of  minimum  disturbance  of  the  plasma  it  would  be 
very  attractive  to  be  able  to  couple  rf   energy  into  a  fusion 
plasma   at  technically  reasonable  frequencies.   The  main  contri- 
bution of  the  present  work  to  the  subject  of  plasma  heating  is  to 
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document  computations  which  support  the  assertion  that  rapid  rf 
energy  absorption  and  eventual  plasma  heating  in  thermonuclear 
devices  may  be  produced  at  frequencies  several  orders  of  magni- 
tude less  than  the  ion  cyclotron  frequency  through  the  excitation 
of  Alfven  waves.   A  secondary  contribution  is  to  provide  a 
comprehensive  review  of  the  theoretical  foundations  of  this 
absorption  and  heating  mechanism. 

The  heating  of  plasmas  by  Alfven  waves  was  actually  proposed 
earlier  by  Furth  (1959) .   This  early  proposal  had  the  disadvan- 
tage of  requiring  that  the  wave  generator  be  imbedded  in  the 
plasma,  a  condition  which  is  clearly  unreasonable  for  thermo- 
nuclear plasmas.   Recently  though,  several  authors  (Grossmann 
and  Tataronis,  1973;  Tataronis  and  Grossmann,  1974;  Tataronis, 
1975;  Kappraff,  Tataronis  and  Grossmann,  1975;  Tataronis  and  Gross- 
mann, 1976;  Tataronis  and  Kappraff,  1977;  Chen  and  Hasegawa,  1974) 
have  considered  an  alternative  approach  to  heat  plasmas  by  means 
of  Alfven  wave  excitation.   The  physics  of  this  new  approach  can 
best  be  understood  by  considering  the  plasma  to  be  a  perfectly 
conducting  magnetohydrodynamic  (MHD)  fluid  with  profiles  in- 
homogeneous  in  the  direction  perpendicular  to  the  magnetic  flux 
surfaces.   As  shown  elsewhere  (Tataronis  and  Grossmann,  1971, 
1972;  Uberoi,  1972),  Alfven  waves  in  the  MHD  model  propagate 
locally  on  magnetic  flux  surfaces,  and  have  a  phase  velocity 
which  varies  in  space  according  to  the  spatial  variations  of  the 
plasma  profiles.   This  local  variation  in  the  Alfven  phase 
velocity  and  the  fact  that  Alfven  wave  propagation  is  one 
dimensional  along  the  magnetic  flux  surfaces  together  imply  the 
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existence  of  a  continuous  spectrum  of  Alfveh  wave  frequencies. 
In  this  continuum,  the  eigenfunctions  of  the  linearized  MHD 
operator  are  not  square  integrable  due  to  the  existence  of 
spatial  singularities,  and  consequently  energy  can  accumulate 
without  bound  at  the  point  of  singularity.   Among  the  many 
effects  that  can  be  attributed  to  this  accumulation  or  equiv- 
alently  absorption  of  rf  energy,  we  mention  here  damping  of 
MHD  surface  waves.   This  damping  in  an  inhomgeneous  plasma 
column  in  MHD  has  been  predicted  by  several  authors  (Grad,  1969; 
Grossmann  and  Tataronis,  1973),  the  damping  rates  calculated 
for  special  configurations  (Grossmann  and  Tataronis,  1973) 
and  experimentally  verified  (Grossmann,  Kaufmann  and  Neuhauser, 
1973) .   The  experimental  verification  of  the  damping  of  MHD 
surface  waves  was  accomplished  in  a  plasma  in  which  one  would 
have  expected  significant  departures  from  the  ideal  MHD  model, 
and  consequently  adds  strong  support  to  the  assertion  that  the 
proposed  absorption  mechanism  is  physically  real  and  potentially 
important.   More  recently,  direct  experimental  observation  of 
plasma  heating  by  Alfven  waves  has  taken  place  in  several 
laboratories.   Specifically,  Golovato,  Shohet  and  Tataronis  (1976) 
have  observed  heating  in  the  Proto-Cleo  Stellarator,  while  Pochelon 
and   Keller   (1977)  have  observed  this  heating  in  a  theta  pinch. 
The  experimental  measurements  by  Pochelon  and  Keller  are  in 
excellent  quantitative  agreement  with  the  theoretical  predictions 
of  Tataronis  and  Grossmann  (1974,  1976).   Other  experiments  have 
been  carried  out  by  Vo,  et  al  (1976)  . 
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The  first  complete  treatment  of  the  theory  of  resonant 
absorption  in  an  MHD  plasma  in  the  continuous  spectrum  of  stable 
Alfven  oscillations  was  given  by  Tataronis  (1975) .   A  special 
example,  i.e.,  incompressible  linearized  MHD  with  spatially 
varying  planar  plasma  profiles,  was  examined  and  the  basic 
underlying  theory  of  the  plasma  dynamic  response  to  a  time 
varying  external  source  oscillating  at  frequencies  corresponding 
to  the  continuum  of  plasma  Alfve*n  frequencies  was  given.   The 
results  of  the  theory  show  that  surprisingly  good  coupling 
between  the  source  and  the  plasma  can  be  obtained.   The  calcula- 
tions presented  in  the  present  paper  represent  an  extension  of 
the  theory  of  Tataronis  (1975)  ,  and  the  main  new  feature  of  the 
present  work  is  the  complete  treatment  of  currently  realistic 
plasma  configurations  in  cylindrical  geometry  in  a  fully  com- 
pressible ideal  MHD  model.   In  particular,  we  consider  reso 
absorption  in  such  currently  interesting  plasmas  as  ST,  ORMAK 
and  Scyllac;  these  plasma  configurations  are  characterized  by  large 
magnetic  field  shear  and  strong  radial  inhomogeneity.   The  calcula- 
tions thus  are  made  for  both  low  and  high-B  configurations.   One 
of  our  goals  is  to  point  out  that  the  shape  of  the  plasma  profiles 
plays  a  major  role  in  the  quantitative  energy  coupling  and 
absorption. 

Another  feature  of  the  present  calculations  is  that,  due  to 
compressibility,  a  second  continuum  arises  in  ideal  MHD,  and 
this  continuum  also  results  in  resonant  absorption.   The  second 
continuum  is  conventionally  referred  to  as  the  cusp  continuum 
and  is  associated  with  the  singular  one-dimensional  nature  of  a 
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part  of  the  slow  magnetoacoustic  wave.   We  show  that  the  cusp 
continuum  indeed  yields  positive  absorption  when  the  external 
exciting  frequency  lies  in  the  required  frequency  range  but  does 
not  represent  a  potentially  interesting  possibility  for  heating. 
The  reasons  are  twofold.   First,  in  all  cases  investigated  here 
the  Alfven  continuum  yields  much  stronger  absorption  and  second, 
the  cusp  continuum  disappears  when  the  guiding  center  model  is 
used  to  describe  more  accurately  the  plasma  (Tataronis  and 
Grossmann,  1976)  . 

Finally,  as  will  become  evident  later,  the  numerical  results 
indicate  that  very  strong  absorption  can  be  found  for  practically 
all  types  of  plasma  profiles  for  exciting  frequencies  lying  in 
the  Alfven  continuum.   These  results  warrant  attention  since 
the  external  apparatus  necessary  to  produce  the  required  energy 
source  is  modest. 

The  remainder  of  the  paper  is  divided  into  eleven  sections. 
Section  2  presents  the  basic  equations  and  mathematical  details 
concerned  with  the  ideal  MHD  model.   In  that  section,  the  basic 
energy  integral  upon  which  the  absorption  theory  is  based  is 
derived.   Section  3  connects  the  classical  folklore  of  resonant 
absorption  in  homogeneous  plasmas  to  that  of  the  present  in- 
homogeneous  case.   Section  4  considers  pertinent  characteristics 
of  the  solutions  to  the  basic  equations.   Such  features  as  the 
singular  nature  of  the  eigenmodes  and  their  representation  are 
discussed.   Section  5  develops  an  asymptotic  form  of  the  energy 
absorption  rate  as  a  function  of  the  solutions  obtained  in 
Section  4.   A  representation  of  the  effective  coil  impedance  is 
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developed  in  Section  6.   The  actual  plasma  profiles  used  are 
displayed  and  discussed  in  Section  7,  and  numerical  techniques 
employed  to  solve  for  the  impedance  for  these  profiles  are 
discussed  in  Section  8.   The  results  of  the  numerical  computa- 
tions are  given  in  Section  9.   Section  10  is  devoted  to  a 
summary  of  the  experimental  work  on  Alfven  wave  heating  that  has 
been  accomplished  to  this  date.   The  conclusions  of  this  paper 
are  given  in  Section  11. 
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2 .   Basic  Equations 

In  this  section  we  present  the  pertinent  mathematical  details 

concerning  the  calculation  of  the  energy  absorption.   It  should 

be  pointed  out  that  the  mathematical  treatment  follows  very  closely 

that  of  Tataronis  (1975)  and  differs  primarily  in  the  use  of 

cylindrical  coordinates   which  allow  realistic  magnetic  field 

shear.   A  further  difference  arises  in  the  consideration  of  fully 

compressible  MHD  modes.   The  plasma-vacuum-source  system  which 

we  consider  is  shown  schematically  in  Figure  1.   A  perfectly 

conducting  diffuse    plasma  extends  to  a  radius   r  where  it  is 

joined  to  a  pure  vacuum  region.   An  ideal  sheet  current  is  placed 

at  radius  r   ,  and  the  entire  system  is  surrounded  by  a  conducting 

wall  at  radius  r  . 
w 

The  Governing  Equations 

The  plasma  is  described  at  this  point  by  generalized  spatially 
varying  profiles  for  the  magnetic  field,  pressure  and  density 
respectively  as : 

B  =  [0,BQ(r) ,Bz(r)]  ,  (1) 

P  =  P(r)  ,  (2) 

P  -  P(r)  ,  (3) 

where   B   and  B   are  the  6  and  z  components  of  the  equilibrium 
magnetic  field.   All  quantities  vary  with  the  radius   r   and 
as  r  ♦  r   ,   P(r)  and  p(r)  vanish  smoothly.   The  equilibrium 
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relationship  of  ideal  MHD  connects  P(r)  and  B(r)  through  the 
following  equation, 

a  r2     Bl 

-5-2-  (P  +  JL.)  +_1=0/  (4) 

dr      2\xQ'         rpQ  v  ; 

2 
where  y.  is  the  vacuum  permeability  and  B   has  been  written 

2     2 

for  Bg  +  B  .   In  the  vacuum  an  externally  supported,  helical 

surface  current  is  placed  at  the  position  r  =  r  and  is 
represented  as, 

Js(r,9,z,t)  =  [Jse(9,z,t)9  +  Jsz(6,z,t)2]  6(r-rc)  ,   (5) 

where  6  (r-r  )  is  the  Dirac  delta  function  and, 
c 

J  fl  =  Jn(t)  cos  v  cos  (m  6  +  k  z)  ,  (6) 

J   =  Jn(t)  sin  v  cos  (m  9  +  k  z)  ,  (7) 

s  z     u  c      c 

tan  v  =  -  m  /k  r  .  (8) 

Here,  v  is  the  pitch  angle  that  the  surface  current  makes  with 

a  plane  perpendicular  to  the  z-axis  as  shown  in  Fig.  1,  and  J. 

and  J   satisfy  the  following  condition, 
z 

mc  j  _  +  k  J   =  0  (9) 

—  us9     c  sz 

r 
c 

We  point  out  that  J   has  the  units  of  current  per  unit  area, 

while  Jg  and  J   are  measured  in  terms  of  current  per  unit  length, 

length  being  measured  in  the  direction  perpendicular  to  the 

current  density  on  the  surface  r  =  r  •   This  surface  current  is 
J  c 

the  source  of  the  external  force  which  acts  on  the  plasma  surface. 
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Practically,  a  surface  current  having  helical  symmetry  can  best 
be  approached  by  line  currents  which  twist  around  the  plasma 
column,  similar  to  that  in  the   Stellarator  configuration. 
However,  due  to  the  discrete  nature  of  such  line  currents, 
harmonics  of  the  specified  m  and  k  values  are  generated  in  the 
device.   Shohet,  Golovato  and  Tataronis  (197  6)  have  shown  that 
these  harmonics  have  no  adverse  effects  on  the  absorption 
mechanism  under  consideration  here.   On  the  contrary,  the  harmonics 
may  lead  to  the  desirable  condition  of  local  energy  absorption 
at  several  points  across  the  plasma  column. 

The  plasma  dynamics  are  described  by  the  linearized  ideal  MHD 
equations  of  motion  which  allow  adiabatic  compressible  displace- 
ments. The  equations  combined  with  Maxwell's  equations  reduce 
to  the  following  set, 

8Y   i  i 

PJT"  -*-  (VxB)xQ  +  -±_  (VxQ)xB  "  Vp  ,  (10) 

ot      uQ  ~       -      MQ      -   - 

3Q 

*£  =  (B-V)v  -  (v-V)B  -  B  V-v         ,         (11) 

9p 

_  =  -  (vV)P  -  TP  V-v  ,         (12) 

where   v,  Q,  p   represent  the  linearized  perturbation  velocity 
vector,  magnetic  field  vector  and  pressure,  respectively.  The 
symbol  T    is  the  ratio  of  specific  heats.   All  other  variables 
represent  previously  defined  equilibrium  quantities. 

In  the  vacuum  region  the  fluctuating  perturbed  electric, 
<5E   ,  and  magnetic,  6B   ,  fields  are  determined  through, 
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6Bv   =   V<j>              ,  (13) 

ft    6?v   =   "    V    X6?v    '  {14) 
where   the  magnetic   potential   $   satisfies, 

V2<j>    =   0    .  (15) 

At  r  =  r   ,  the  plasma  and  vacuum  fields  are  connected  by 
the  usual  continuity  conditions  of  ideal  MHD  which,  in  linearized 
form,  reduce  to, 


n  x  6Ev  =  (n-v)Bv  ,  (16) 


-  TPV-5  +  -£-  B-[Q  +  (£-V)B]  =  -L.  B  •  [SB   +  (£-V)B  ]  ,    (17) 


where  C(r,t)  is  the  Lagrangian  displacement  vector  related 
to  v   through, 


Ur,t)  = 


dt*  v(r,f)  .  (18) 


Equation  (16)  requires  Q   to  be   parallel  to  the  fluctuating 
plasma-vacuum  interface,  and  eq.  (17)  requires  the  jump  in 
total  pressure  across  the  fluctuating  interface  to  vanish.  The 
vector  n  is  the  unit  vector  normal  to  the  interface.  The 
subscript  (  )   denotes  vacuum  quantities.   Finally,  the  following 
jump  conditions  across  the  surface  current  at  r  =  r  must  be 


satisfied, 


n  x  [6Bv]  =  yQJs  ,  (19) 

n  x  [6Ev]  =  0     ,  (20) 
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as  well  as  the  condition  that  the  perturbed  magnetic  field  be 


parallel  to  the  metallic  wall  at  r  =  r   ,  i.e., 

n*6B   ■  0  . 


(21) 


Expressed  in  terms  of  cylindrical  variables,  these  constraints 
reduce  to  the  following  system: 

r  -»■  0  : 

(i)       all  dependent  variables  remain  bounded; 


(ii) 


r  =  r   : 
P 

U0  ~v   -v      yQ    r 


36B. 
r    3t 


vr 


(22) 


(23) 


r  =  r 


(iii)   • 


[6Bvz]  =  "  Vse  ' 


[6B   ]  =  0 
vr 


(24) 


(25! 


r  =  r   : 
w 


(iv) 


6B    =  0 
vr 


(26) 


The  brackets  in  the  above  equations  denote  the  jump  in  the 
enclosed  quantities  in  the  direction  of  increasing  r. 


•11- 


In  a  manner  similar  to  that  of  Tataronis  (1975) ,  we  make  use 
of  the  displacement  vector  £  to  derive  a  relationship  which 
expresses  the  rate  of  change  of  the  plasma  energy  in  terms  of 
the  fields  induced  by  the  surface  current.   The  procedure  for 
this  derivation  involves  integrating  eqs.  (11)  and  (12)  with 
respect  to  t,  with  the  initial  conditions  set  equal  to  zero, 
substituting  the  resulting  expressions  for  Q  and  p  as  functions 
of  £  into  eq.  (10),  scalar  multiplying  eq.  (10)  by  v  and  then 
integrating  over  a  plasma  region  bounded  by  a  surface  Z.   The 
following  result  is  obtained, 


dW 

dt 


=  -  o  do-v  p*   ,  (27) 

J   ~  ~ 

Z 


where  p*  has  been  written  in  place  of  p  +  B*Q/u0,  the  total 
perturbed  plasma  pressure,  and  W   represents  the  sum  of  the 
plasma  kinetic  and  potential  energies  induced  by  the  linearized 
fields ,  i.e., 


«p  =  I 


dx[— q2  +  j-  (CxQ)  +  (v-£K-vp+rp(v-c)2] 


2    1  f 
dx  pv  +  - 

2  J   kW0  - 

V  V 

P  P  (28) 


Here,  V   is  the  plasma  volume  enclosed  by  Z.   In  what  follows, 

V  will  be  defined  by  r  =  rn  ,  0  <  8  <  2tt  and  0  <  z  <  L   (E  2-n/k   ) 
p  J      P      -    -  -    -   z        c 

with  periodic  boundary  conditions  imposed  on  the  perturbed 
variables  at  z  =  0  and  z  =  L  as  well  as  at  9  =  0  and  6  =  2tt. 
Equation  (27)  states  that  the  rate  at  which  the  plasma  energy 
increases  with  time  is  equal  to  the  energy  flux  density,  vp*  , 


-12- 


directed  into  the  volume  integrated  over  the  enclosing  surface 
Z .   If  E  is  identical  with  a  plasma-vacuum  interface  where  the 
boundary  conditions  given  by  eqs .  (16) -(17)  apply,  then  it  is 
a  straightforward  matter  to  express  the  right-hand  side  of 
eq.  (27)  in  terms  of  the  component  of  the  vacuum  Poynting 
vector  normal  to  I,  which  in  turn  can  be  related  to  the  time 
rate  of  change  of  the  energy  associated  with  the  perturbed 
vacuum  magnetic  field  and  to  the  power  delivered  by  the  surface 
current.   For  the  geometry  in  Figure  1   and  under  the  assumption 
that  B„ (r)  varies  continuously  across  the  equilibrium  plasma- 
vacuum  interface,  the  following  representations  of  dW  /dt  are 
readily  obtained: 


dW 


2tt 


dt 


_E  = 


"  frPd9f 

0        0 


dz  v  (r  ,t)p  (r  ,t)  , 
r  p       p 


(29) 


2tt 


L 


rp  de 


dz  v  (r,t) 
r  P 


B  -6B  (r  .t)  u  1 
~v  -v  p     o 


(30) 


2  it 


rp  d6  I   dz  Sr(rp,t)  , 


(31) 


dW 


2tt 


v 


dt 


L 


r   de 
c 


dz  6Ey(rc,t) -Js  , 


(32) 


where   S   is  the  radial  component  of  the  vacuum  Poynting  vector, 


6E   x  <5B   ,  and, 


f  6B* 

wv  -  J  dTv  rir  ' 


(33) 
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W  being  the  energy  of  the  perturbed  vacuum  magnetic  field  and 
dx   the  differential  volume  element  taken  in  the  vacuum. 

Equations  (10)  -  (33)  form  the  basis  of  our  investigation. 
In  order  to  describe  the  absorption  mechanism  which  is  deduced 
from  these  equations,  it  is  useful  to  make  use  of  certain  local 
properties  that  characterize  the  linear  MHD  waves  in  infinite 
uniform  plasmas.   A  discussion  of  these  waves  is  therefore 
given  in  the  next  section  where  it  is  shown  that  two  of  the 
three  MHD  modes  imply  local  motion  of  the  plasma.   It  is  this 
local  motion  that  is  the  source  of  the  continuous  spectrum  of 
the  MHD  differential  operator  and  of  the  resulting  energy 
absorption.   In  later  sections  we  shall  solve  eqs.  (10) -(12) 
for  v  (r,8,z,t)  and  p  (r/6,z,t)  in  terms  of  the  specified 
surface  current  and  use  these  solutions  to  obtain  an  explicit 
representation  of  dW  /dt  through   eq.  (29) .   We  make  use  of  the 
resulting  expression  to  derive  the  experimentally  relevant  coil 
impedance  and  then  evaluate  this  impedance  for  Tokamak-  and 
Scyllac-like  plasmas. 
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3 •   Resonances  in  the  Infinite  Homogeneous  Plasma 

In  this  section  we  will  briefly  review  MHD  wave  propagation 
in  infinite,  homogeneous  plasmas.   Specifically,  we  will 
connect  the  classical  concepts  of  resonance  and  cut-off  frequen- 
cies in  the  dispersion  characteristics  of  the  MHD  waves  to  the 
subject  of  the  present  paper. 

We  assume  that  the  equilibrium  magnetic  field  is  unidirec- 
tional  throughout  the  plasma  and  that  the  equilibrium  fields 
do  not  vary  in  space.   Assuming  propagation  in  the  form 
exp  i  (cot  -  k-r),   one  can  readily  derive  from  eqs.  (10)-(12)  the 
following  dispersion  relation  for  the  linearized  MHD  waves: 


-2.22    2„.  ,  4   ,2     2,,  2  2     2  2,  4    2on      .      ,,., 

(id  -k  v  cos  8)  [co  -  (v   +  v  )  k  lo   +  v  v  k  cos  8]   =   0  ,    (34) 

a  3,       S  cl  S 


where   k  represents  the  magnitude  of  the  wave  vector  k,  8  is 
the  angle  between  k  and  B,  and  v   (E   B//pu~7)  and  v   (=  /fP/p) 

mm  -~  A  US 

are  respectively  the  Alfv£n  and  sound  speeds.  Equation  (34) 
yields  three  modes  of  propagation:  the  Alfv£n  wave, 


2  2    2 
cj   -  k  v  cos  8  =  0  ,  (35) 

& 


and  the  slow  and  fast  waves, 

4          i    2            2.  .2    2  2    2.4         2a         n                              ,,-, 

to      -(v      +v)kio  +vvk   cos    0    =    0.                         (36) 

cl       S  3  S 

An  illuminating  representation  of  the  solution  of  eqs.  (35)  and 
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(36)  is  found  in  the  phase  velocity  diagram  (the  Friedrichs 
diagram)  shown  in  Figure  2,  where  the  loci  of  tu/k  is  traced 
as  a  function  of  0.  The  envelopes  have  been  drawn  assuming 
v   <  v-  .   For  v   <  v   ,  the  envelope  of  the  Alfve'n  is  tangent 

So.  S3. 

to  the  fast  wave  envelope  at  6  =  0.   One  sees  in  Fig.  2  that 
the  phase  velocity  of  the  fast  wave  is  more  or  less  isotropic 
in  9  but  that  of  the  slow  and  Alfven  waves  is  highly  anisotropic. 
Of  particular  importance  is  the  fact  that  as  Q   approaches  90°, 
corresponding  to  perpendicular  propagation,  the  magnitude  of  the 
phase  velocity  of  the  slow  and  Alfve'n  waves  approaches  zero. 
The  implication  of  this  property  is  that  localized  excitation 
of  the  fluid  gives  rise  to  perturbations  which  propagate  local- 
ly, i.e.  free  of  dispersion  in  the  direction  perpendicular  to 
B,  about  the  magnetic  line  of  force  on  which  the  initial  data 
was  specified.   This  fact  is  borne  out  in  Fig.  3  where  the  wave 
fronts  formed  by  the  three  MHD  waves  excited  by  a  source  localized 
at  the  origin  of  an  x-z  Cartesian  coordinate  system  are  plotted 
at  some  instant  of  time.   The  wave  front  of  the  fast  wave  is 
isotropic  as  one  would  expect,  while  those  of  the  slow  and 
Alfven  waves  are  localized  to  the  interior  of  a  cone.   The  Alfven 
wave  is  actually  degenerate  in  the  sense  that  it  has  no  dispersion. 
As  indicated  in  Fig.  3,  a  point  disturbance  in  the  Alfven  mode 
remains  a  point  disturbance  localized  to  the  field  line  on  which 
the  excitation  was  made.   Again  the  Alfven  wave  in  Fig.  3 

represents  the  case  v   <  v  .   For  v   >  v  ,  the  Alfven  mode  is 

s    a        s     a 
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changed  only  by  moving  the  point  representing  the  wave  front 
from  the  slow  wave  envelope  to  the  fast  wave  envelope.   The 
slow  wave  also  reveals  localization  in  the  form  of  the  vertex 
of  the  cusp  trailing  to  the  left  of  the  wave  front.   Considerable 
effort  has  been  made  in  the  past  to  investigate  the  mathematical 
and  physical  consequences  of  the  "cusp"  motion  (Cumberbatch, 
1962).   The  main  results  are,  (1)  the  cusp  origin  is  a  point 
which  moves  along  the  magnetic  field  line,  (2)  the  cusp  motion 
is  a  one-dimensional   motion,  and  (3)  the  speed  at  which  the  cusp 
moves  is  given  by, 

2         1       rP    b2/^0 

VC     p r-2-  (37) 

c      p  rp  +  bVu0 

The  localization  of  the  Alfven  and  cusp  propagation  is  the  origin 

of  the  MHD  resonances  that  we  shall  refer  to  later  as  the  cusp 

and  Alfven  wave  resonances.   In  order  to  bring  this  resonance 

phenomenon  more  in  line  with  the  terminology  used  in  plasma  wave 

theory,  we  show  in  Fig.  4  the  dispersion  diagram  for  the  three  MHD 

waves.   There,  the  frequency  w  is  plotted  as  a  function  of  kj_  , 

the  component  of  the  wave  vector  perpendicular  to  B,  with  k   , 

the  parallel  component,  held  fixed.   In  wave  theory  a  resonance 

is  said  to  occur  where  the  wave  vector  becomes  infinite.   In  Fig. 

4  two  such  resonances  are  found.   One  occurs  at  u  =  m,  =  k,v  ; 

a    ii  a    a 

will  be  referred  to  from  this  point  on  as  the  Alfven  frequer. 
Another  such  resonance  occurs  at  u  ■  u.    k  v  ;  w  will  similarly 
be  referred  to  as  the  cusp  frequency.   One  also  notices  in  Fig.  4 
that  with  k   held  fixed,  the  three  MHD  waves  have  cut  off 
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frequencies,  i.e.  frequencies  at  which  the  wave  number, k± 
in  this  case,  vanishes. 

The  importance  of  resonance  and  cutoff   frequencies  arises 
in  connection  with  the  prediction  of  energy  absorption  and 
reflection  in  spatially  inhomogeneous  plasmas.  Specifically, 
if  one  considers  an  inhomogeneous  equilibrium,  a  differential 
equation  replaces  the  wave  dispersion  relationship.  The  solu- 
tions of  the  differential  equation   are  characterized  by 
singularities  which  occur  in  space  at  the  point  where  the 
applied  frequency  equals  the  local  value  of  the  now  spatially 
dependent  resonance  frequency.   Thus,  one  finds  a  correspon- 
dence between  the  resonances  of  an  infinite  uniform  plasma  and 
the  spatial  singularities  of  certain  macroscopic  fluid  variables 
in  an  inhomogeneous  plasma.   Depending  on  the  nature  of  the 
singularities  and  the  dependence  of  the  wave  energy  on  the 
perturbed  fields,  one  can  establish  in  certain  cases  an  associa- 
tion between  these  singularities  and  unbounded  energy  at  the 
singular  point.   Similarly,  one  can  demonstrate   a  connection 
between  cutoff   frequencies  in  the  dispersion  characteristics 
of  the  homogeneous  medium  and  energy  reflection  in  the  inhomo- 
geneous medium.   Based  on  these  qualitative  ideas,  one  can  now 
make  the  following  statements  in  regard  to  energy  absorption  in 
an  ideal  MHD  fluid,  namely  that  energy  absorption  may  be  possi- 
ble at  the  local  Alfv£n  wave  and  cusp  resonances.   In  the 
following  sections  we  shall  demonstrate  the  validity  of  this 
assertion.   We  mention  here  that  other  rf  plasma  heating  schemes, 
such  as  ion  and  electron  cyclotron  heating  and  lower  hybrid 
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resonance  heating,  are  based  on  wave  resonances  of  this  type. 

The  above  discussion  has  been  concerned  with  an  infinite  homo- 
geneous plasma;  the  main  point  of   the  present  paper  is  that 
it  is  not  only  useful  but  absolutely  necessary  to  consider  the 
inhomogeneous  nature  of  real  plasmas  in  order  tc  treat  correctly 
the  absorption  of  rf  energy  in  a  plasma  by  means  of  such  reson- 
ance phenomena.   In  fact  the  results  which  we  will  later  present 
show  that  the  questions  of  existence  of  absorption,  the  magni- 
tude and  the  position  in  the  plasma  of  maximum  absorption  depend 
critically  on  the  detailed  shape  of  the  inhomogeneous  plasma 
profiles  such  as  those  of  the  pressure,  density  and  magnetic  field, 
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4 .   Solution  of  the  Equations  of  Motion 

In  this  section  we  begin  our  analysis  of  the  linearized  MHD 
equations  in  cylindrical  geometry.   In  particular  we  shall  solve 
eqs.  (10)  -  (12)  in  the  presence  of  a  time  varying  forcing  term 
and  then  apply  the  results  to  compute  the  rate  at  which  energy 
is  transferred  to  the  plasma  as  a  result  of  the  action  of  the 
forcing  term.   Moreover  we  will  discuss  certain  pertinent 
details  of  the  spectrum  of  these  equations  and  discuss  the 
nature  of  solutions  in  the  vicinity  of  singularities.  This 
information  will  then  be  useful  in  describing  the  energy  absorp- 
tion processes  in  the  inhomogeneous  plasma. 

For  the  present  calculation,  there  is  a  natural  choice  of 
macroscopic  variables  to  use  in  solving  the  MHD  equations . 
These  variables  are  the  total  perturbed  pressure  p  ,  where 
P  =  P  +  Q*B/yn   ,  and  the  radial  perturbation  displacement 
vector  £  .   These  variables  occur  explicitly  in  the  linearized 
expressions  for  the  boundary  conditions  at  the  vacuum  interface. 
Further,  the  governing  differential  equations  for  these  two 
variables  greatly  facilitates  an  investigation  of  the  linearized 
MHD  spectrum  as  shown  by  Grad  (1974) . 

We  first  perform  a  Laplace  transformation  with 
respect  to  time  and  a  Fourier  transform  with  respect  to  6  and  z  by 
writing, 


5r(r,6,z,t)  = 


ao 

dw  r 

2tt 

—  ao 


^   V    exp  i(ut-me-kz)  E,    (r,m,k,u))  , 
m=-°° 

Im  (to)  <  0,     (38) 
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where 


£r  (r,m,k,uj)  = 


dt 


f 


2tt 


dz 


■^   exp[-i  (d)t-m8-kz)Cr  (r,e  ,z,t)  , 


(39) 


L   being  a  contour  in  the  complex  (j  plane  below  all  singularities 
of  the  integrand.   The  other  dependent  variables  are  to  be 
treated  in  an  analogous  manner.   After  application  of  eq.  (39), 
a  straightforward  manipulation  of  eqs .  (10)  -  (12)  yields  the 
following  two  differential  equations  for  p   and  £   : 


dp      ^ 

gjr-  =  bQp   +  blrlr    , 


(40) 


dr^r 

3r—  =b2P   "  V^r 


(41) 


where  we  have  set  all  initial  conditions  equal  to  zero. 
Equations  (40)  and  (41)  can  be  combined  to  yield  a  single 
second  order  equation  for  £   : 

,2 


d_ 
dr 


1   dr^r 


b2   dr 


bi  +  b^  "  dT 


ri      =   0  . 
r 


(42) 


The   coefficients      b_,    b.    and   b-    are   defined   as    follows.    Let, 


2        F* 
A    =    -    p    0)       +    —      , 

Uo 


c  ■  -   (rp    + 


B            ^  ^2 

2  F 

)u  p  +  rp  —    , 

yo  uo     ' 


(43) 


(44) 


D   -    ui4p2    +       (51y   +    k2)C    , 

r 


(45) 
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2     2     2 

where  as  previously  defined  B   =  Bq    +   B  ,  and, 


F  =  -  BA  +  k  B 
r  H       z 


(46) 


The  coefficients  b_  ,  b,  and  b„  then  can  be  written  as, 


1_ 

b. 


AC 

rD 


(47) 


b0  =  "  I 


f,    V 

' 

I  rv 

m  f  +  r    u4p2 

F  F  +  Be  ^c^ 


k     i  fa  o  Bq   d  (Bel 

b,  =  +  —  A-2  r  -7—  —  - 

1     r  |_    ryn    dr  [r   t 


Ba^2  F2 


A 


ru, 


(48) 


,  Be]2  *4p2    i^ 

2    -r—     „   B 


rUQJ    AC   D9 


(49) 


Equations  (40) -(42)  have  been  used  by  many  authors  (Grad, 
1973;  Goedbloed  and  Sakanaka,  1974)  to  discuss  the  MHD  spectrum 
and  the  normal  mode  approach  to  plasma  stability  of  the  screw 
pinch.   The  absorption  process  under  discussion  here  depends  on 
the  singular  nature  of  the  above  equations.   Specifically,  these 
equations  are  singular  at  points  r  =  rfi ( w)  where  the  coefficients 
bQ,  b,  or  b_  possess  singularities,  i.e.,  at  the  zeros  of  A(r,w) 
and  C(r,u>).   The  factors  A  and  C  set  equal  to  zero  correspond 
respectively  to  the  local  Alfven  and  cusp  dispersion  relations 
of  the  infinite  homogeneous  plasma  considered  in  Section  3,  but 
written  here  in  terms  of  cylindrical  coordinates.   The  zeros 
of  A  satisfy, 


2     2.  ,    1    F*(r) 

a)  =  ul  (r)  =  —  ) — (-    , 

^      M0    P(r) 


(50) 


while  those  of  C  are  obtained  from, 
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to2  =  w2(r)  ,  _1 rP(r)F2(r)  (51) 

V(r>   rP(r)+B2(r)/u0 


where  co   and  w   are  respectively  the  local  frequencies  of  the 
Alfven  wave  and  cusp  resonances  defined  in  Section  3.   Thus 
about  any  r  =  rQ(oj)  satisfying  eqs.  (50)  or  (51),  the  solutions 
of  eqs.  (40)  and  (41)  may  be  singular.   By  expanding  all 
coefficients  in  a  power  series  about  the  point  r-(u)),  one  can 
prove  the  existence  of  a  singular  solution.   One  finds  in  fact 
the  following  two  linearly  independent  solutions  to  eq.  (42), 


f   =  g(r,u))  (52) 

f   =  g(r,u>)  In    (r.  -  r)  +  (r  -  r  )  h(r,uj)  ,       (53) 


where  g(r,u)  and  h(r,w)  are  analytic  in  some  domain  about  r_, 

the  extent  of  the  domain  being  determined  by  the  closest  other  zero 

of  A(r,a>)  or  C(r,io)  to  rQ.   The  general  solution  for  £   can 

thus  be  written  as  follows, 

rC   =  C  f   +  C  f    ,  (54) 

r    a  a    s  s 

where  C  (oj,m,k)  and  C  (io,m,k)  are  functions  independent  of  r 

cl  S 

that  are  determined  by  invoking  boundary  conditions  at  r  =  0 
and  r  =  r  .   Hence,  for  some  given  u>,  the  variable  £   possesses 
in  general  a  logarithmic  singularity.   Actually  the  singularity 
manifests  itself  in  a  manner  stronger  than  a  logarithm.   This 
can  be  seen  in  the  following  manner.   Let  us  first  solve  for  £9 
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and  £  ,  the  azimuthal  and  longitudinal  components  of  t   respec- 
tively  in  terms  of  r£   and  dr£  /dr  .   A  lengthy  calculation 
yields: 


2   Bz      m    drCr 

-  irD£e  =  (-  pu,2E  ^  +  TP  S  A)  -g-E 

+  2  -i-  (-poj2B   +  TPkF)kr£    ,  (55) 

->   BQ  dr£ 
-irDCz   =  (-  pu/E  —  +  TPkA)  -^ 

B6      2         m 

"  2  ^  (-pw  BQ  +  TP  g  F)kr£r   ,  (56) 


where/ 

E  =  F  Bz  ~  kB9  *  (57) 

The  apparent  singularity  at  the  zeros  of  D  is  misleading,  for 

it  can  be  demonstrated  in  terms  of  the  solutions  for  ?   that 

sr 

£g  and  i      are  analytic  at  points  where  D  vanishes.   By  substi- 
tuting eq.  (54)  for  ?   in  eqs .  (55)  and  (56),  one  can  readily 
show  that  about  rf)(o)), 

^e  and  K   *  F^T  •  (58) 

These  variables  possess  stronger  singularities  than  those  of 
r£.   in  the  sense  that  t,   and  £   are  not  square  integrable  whereas 
r£   is  square  integrable.   As  pointed  out  in  the  introduction  and 
as  will  be  expanded  upon  later,  the  non-square  integral  singu- 
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larities  of  £„  and  E,      form  the  basis  of  the  energy  absorption 
process  which  is  the  subject  of  the  present  paper. 

In  order  to  proceed  with  a  calculation  of  the  time  behavior 
of  £  ,  it  is  necessary  to  consider  the  coupling  of  the  sheet 
current  to  the  plasma  and  determine  the  singularities  of 
£  (r,8,z,a>)  in  the  complex  cj  plane,  where  £  (r,0,z,cj)  is  the 
inverse  Fourier  transform  of  %    (r,m,k,u).   The  inverse  Laplace 
integral  can  then  be  performed  on  a  contour  that  lies  below  the 
singular  points.   We  now  briefly  describe  our  construction  of 
E,    (r,e,z,w)  and  then  consider  in  some  detail  the  distribution  of 
the  singularities  in  the  complex  w  plane.   The  inversion  of  the 
Laplace  integral  will  be  considered  at  the  end  of  this  section. 

Let  us  introduce  a  function  'Mr,m,k,a))  which  is  a  solution 
of  eq .  (42),  satisfies  the  boundary  conditions  imposed  on  ?   at 
r  =  0,  namely,  cfj  remain  bounded  in  the  limit  r  ■+   0 ,  and  is 
normalized  by  specifying  di^/dr  at  r  =  0.   It  is  clear  that  these 
conditions  uniquely  define  iji  which  can  then  be  generated  by 
solving  eq .  (42)  as  an  initial  value  problem.   Except  for  special 
equilibrium  profiles,  ij;  must  be  obtained  numerically,  We  now 
write , 

^(r,m,k,io)  =  R(.j,m,k)  ty  (r,m,k,k>)   ,  (59) 

and  determine  R  by  invoking  the  boundary  conditions  at  the  plasma- 
vacuum  interface,  namely,  eqs.  (22)  and  (23).   Assuming  continuity 

of  B0  and  B   at  r  =  r  ,  and  expressing  6B   in  terms  of  the 
o       z         p 

magnetic  potential  i>    through  eq.  (13),  these  conditions  can  be 
expressed  in  the  following  form: 
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at  r   =    r 

E 


rp      d$/dr  b2      ^   dr  0 


and, 


r        ~ 

= £  li 

iFip   dr 


(61) 


where  Fourier  and  Laplace  transformations  have  been  applied, 
and  p*  has  been  eliminated  with  eq.  (41) .   We  now  express  the 
magnetic  potential  <J>  in  terms  of  the  6-component  of  the  external 
surface  current  by  solving  Laplace's  equation  in  the  vacuum 
domain  subject  to  the  jump  and  boundary  conditions  given  by 
eqs .  (24)-(26).   We  give  the  details  of  this  computation  in 
Appendix  A  and  write  the  result  here  as  follows: 

$  =  -  sgn  (k)  i  uQ  rc  JsQ  [Im(|k|r)  -  Xl  Kta(|kjr)]x2  .       (62) 

where, 


Xl 


I  (T  )  -  |k|q  I'(T) 

-HI— E !_! S_P_    ,  (63) 

K  (x  )  -  Iklq  K' (T  ) 
m  p     '  '  ^  m  p 


wjWc)  -  i;(Tc)isk(Tw) 

Im(Tw)  -  *1  Km(V 


(64) 


Here,  I   and  K   are  respectively  the  modified  Bessel  functions 

mm        c  J 

of  the  first  and  second  kind,  the  symbol  (  ) '  designates 

differentiation  with  respect  to  the  argument  of  the  enclosed 

function,  t  ,  x   and  x   have  been  written  in  place  of  (Iklr  )  , 
p    c      w  P 

(|k|r  )  and  ( | k | r.  )  respectively,  k  has  been  taken  as  real,  and 
q,  which  is  related  to  the  surface  impedance  at  the  plasma  inter- 
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face,    is   defined   as    follows, 

q      =      -^ —     at      r   =   r        .  (65) 

d<i>/dr  p 

Finally,  m  is  to  be  interpreted  as  an  absolute  value,  i.e.  as  |m|. 
The  expression  that  is  to  be  substituted  in  eq.  (62)  for  J 

S  o 

is  the  Laplace  and  Fourier  transformations  of  eq.  (6),  i.e., 

Js9  =  tt  J0(U)  cos  v  [6     6(k-kc)  ♦  6      6(k+kc)]   ,    (66) 

c  c 

where  6      and  5(k±k  )  are  respectively  the  Kroniker  and 
m ,  ±m  z 

c 

Dirac  delta  functions.   In  the  analysis  below,  we  shall  assume 
a  sinusoidal  time  dependence  for  J  „  at  the  frequency  wn ,  i.e., 
JQ(t)  =  J   sin  a)  t  from  which  we  obtain, 

Jo(U))  =  Jo  -rr   •  (67) 

aJ0-W 

At  this  point,  we  have  sufficient  information  to  construct 
the  pertinent  dependent  variables  in  terms  of  the  known  function 
^(r,9,k,w).   Equations  (60)  and  (65)  are  combined  to  obtain  q 
which  in  turn  is  substituted  into  eq.  (63)  and  (64)  to  form  Xi/ 
X2  and,  with  eq .  (62),   $.   $  is  then  substituted  in  eq.  (61)  to 
form  R  and  the  result  is  combined  with  eq.  (59)  yielding, 


*«r   =   sgn(k)porcJg9(w)    ^  -$lE$-      ,  (68) 

where 

I' (t    )K' (t     )-I' (t    )K' (t     ) 
.,._  mwmcmcmw 

I'  (T„)K     (Tn)-I     (T)K'  (Tu)-|k|q(u)  [I'll     )K'  (T    )-I'  (Tn)K'  (X     )] 
mwmp        mpmw1  mwmp        mpmw 

(69) 
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The  variables  p*,  5Q  and  £   follow  by  substituting  eq.  (69)  in 
eqs.  (41) ,  (55)  and  (56)  respectively  for  r%    . 

To  obtain  £   in  the  time  domain,  we  first  perform  an  inverse 
Fourier  transformation  of  eq.  (68)  with  respect  to  m  and  k, 
yielding, 

rl    (r,6,z,o>)  =  ynr  J_(u>)  *{"\    f  lr,ai)  ,  cos  (m  9+k  z)  ,   (70) 
r  0  c  0    F(r  )  ij)(r  ,u)      c   c 

where  eq.  (66)  was  substituted  for  J  Q.   In  eq.  (70) ,  k  =  k  , 

so  c 

m  =  m  and  F(r)  =  (m  /r  )BQ(r)  +  k  B  (r) .   The  time  dependence 

C  C   p    o  C  Z 

follows  by  carrying  out  the  inverse  Laplace  integral  of  eq.  (70) 
on  a  contour  lying  below  the  singularities  of  g    (r,6,z,co)  in  the 
complex  co  plane.   The  singularities  are  of  two  types:  poles 
arising  from  the  zeros  of  the  denominator  of  A(co)  and  from  J0(to), 
and  branch  points,  together  with  the  associated  branch  cuts, 
arising  from  the  logarithmic  singularities  of  i|»(r,io)  about  the 
zeros  of  A(r,co)  and  C(r,co).   Details  of  the  branch  point  singu- 
larities can  be  obtained  by  making  use  of  the  technique  employed 
by  Tataronis  (1975)  whereby  in  some  domain,   r.<r<r2,   i|>(r,co)  is 
written  in  the  form, 

*(r,w)  =  Q(co)  [p(co)iJ/  (r,co)  +ipc(r,co)]   ,         (71) 

where  Q(co)  and,p(co)  are  functions  independent  of  r,  and  ^  (r,co) 

a 

and  \l>    (r,co)  are  two  linearly  independent  solutions  of  eq.  (42)  . 

s 

Recalling  that  ip  (r , co)  is  obtained  by  solving  eq.  (42)  as  an 
initial  value  problem  from  r  =  0,  we  assume  that  ^  (r , oo)  has  been 
determined  up  to  a  point  r  =  r..   ip  ( r , oj)  and  dij>  (r  ,co)/dr  at  r  =  r, 
then  serve  as  initial  data  for  the  determination  of  ip(r,co)  in  a 
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domain  r.  <  r  <  r_  where  we  assume  that  either  A(r,<D)  or  C(r,aj) 
vanish  at  r  =  r~  (oo)  for  some  real  id.   We  then  define  ty    (r,uj)  as 
that  solution  of  eq.  (42)  which  is  analytic  at  r  =  r  (to)  and 
if»  (r,w)  as  the  singular  solution  which  we  can  write  as, 


<Ps(r,uj)    1^(^,10) 


^a(r,o))    ^a(r1,co) 


rlDirSuO       dr_!_  _   (?2) 


A(r',u,)C(r',U)  ^2 

r  ra 


Equation  (72)  gives  i/>  (r,to)  in  a  form  whereby  the  relationship 

between  the  zeros  of  A(r',u>)  and  C(r',tjj)  i.e.  r n(u>)  and  the  branch 

point  singularities  of  \p    (r,tu)  in  the  complex  id  plane  for  a  given 

r  is  manifestly  evident.   If  rQ(to)  lies  on  the  real  r'  axis  in 

the  domain  r,  <  r  (to)  <  r,  then  the  Cauchy  integral  on  the  right 

hand  side  of  eq.  (72)  is  singular.   Changing  id  such  that  r  (id) 

moves  off  the  r'  axis  into  the  complex  plane  then  displaces  the 

singularity  of  the  integrand  and  permits  integration  along  the 

real  r'  axis  between  r.  and  r.   Therefore,  contours  in  the  complex 

to  plane  on  which  rn(w)  is  real  and  satisfies  r.  <  r.djj)  <  r  are 

singular  lines  of  the  Cauchy  integral  and  hence  of  ty    (r,co).   Since 

the  coefficients  A(r,w)  and  C(r,uj)  are  real  if  co  and  r  are  real, 

rn(w)  can  be  real  only  if  u)  is  real,  implying  that  the  singular 

lines  of  I/*  (r,oj)  lie  on  the  real  u>  axis.   Typically,  for  real  to, 

r  ('a))  has  the  form  shown  in  Fig.  5  where  the  curve  in  the  domain 

0  <  id  <  (io_)     is  obtained  from  C(r,to)  ■  0,  while  for  id  <  (to.)    , 
C  max  A  mj.n 


we  have  A(r,D)  =  0.   For  (id  )     <  id  <  (<d  )    ,  neither  C(r,j)  nor 

<~  max         A  min 

A(r,<o)  have  real  zeros  for  <d  real.   If  C(r,io)  and  A(r,j)  are 
analytic  in  some  domain  about  the  real  r  axis,  then  the  singular 
lines  of  the  Cauchy  integral  are  branch  cuts.   This  can  be  shown 
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by  assuming  cj  to  be  complex  in  eq.  (72),  w  =  u>   +  i  w.,  and  letting 

<jj .   approach  zero.   One  then  finds  that  if  A(r',<jj  )  =  0  for  some 

real  r'  =  r  ,  between  r..  and  r,  ip  (r,u))  is  defined  on  the  real  w 
axis  as  follows: 


*s(r,(D)    <l>s(rltu) 


t|j^  (r  ,o))    ^  (r,  , oj) 

a.  a    X 


=  P 


r'D(r'  ,u)) 


dr1 


A(r',«i)C(r',»)  #2(rI#Bj 


,  .  D(r.,u))     r. 
+  ie7r  s?"^   ,0         ° 


P<V  C(r0'w)  *2(rnfU) 

3.    U 


dw. 


dr 


-1 


r=r, 


(73) 


where, 


e  =  ■ 


+  1  ,  if  oj .  -+  0  from  negative  values, 


-  1  ,  if  co.  ■+■  0  from  positive  values, 


(74) 


and   P  designates  that  the  principle  value  of  the  integral  is  to 

be  taken,   Therefore,  as  co  cross  the  real  axis  between  u  (r.) 

and  a)  (r)  ,  ij>  (r,cu)  encounters  a  discontinuous  jump  in  value,  imply- 

A  S 

ing  the  presence  of  a  branch  cut  with  branch  points  at  oj  =  ^(r.^) 
and  id  =  oj  (r)  .   Similarly,  if  C(r',m  )  vanishes  at  r1  =  rQ/  one 
finds  the  following  expression  for  i>    (r,co)  on  the  real  co  axis, 


\p s  (r,(jj)         \\>s  (r^,oj) 


ip  a  ( r-  ,oj  )         ^    (r^,oj) 


=  P 


r'D(r'  ,oi) 


dr1 


A<r',u)C<r',«o)    #2(rlfM) 


+   ieir 


sgnU) 


D(ro,w) 


A(r„,u>)    ~2~ 


[rP(r0)+B"(r0)/y0]p(r0)  0'"'    i|»a<r0,u>> 


da)' 


dr 


(75) 


-1 


r=r. 


In  deriving  eqs .  (73)  and  (75),  it  has  been  implicitly  assumed 
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2  2 

that  dw  /dr  and  dco  /dr  do  not  vanish  at  r  =  rQ.   If  the  derivatives 

2  2 

of  wa(r)  and  oj  (r)  do  in  fact  vanish  at  r  =  rn,  one  would  then 

find  that  the  Cauchy  integral  has  a  branch  point  at  u  =  wA(rQ) 

(Tataronis,  1975)  .   In  summary,  we  show  in  Fig.  6  the  expected 

branch  point  and  branch  line  structure  of  ip  (r , oj)  in  the  complex 

2  2 

oj  plane  for  a  given  spatial  distribution  of  ojA(r')  or  oj  (r'). 

2 
Designating  the  resonant  frequencies  as  co.(r')  with  j  =  A  for 

the  Alfven  frequencies  and  C  for  the  cusp  frequency,  we  assume 

2  2 

in  Fig.  6a  a  form  of  oj.(r')  which  has  doj./dr'  =  0  at  r1  =0,  r. 

and  r_ .   As  indicated  in  Fig.  6b,  branch  points  are  therefore 

expected  at  u  =  oj.(O),  ±  oj.(r.),  ±  oj  .  (r_)  and  at  oj  =  ±  oj  .  (r)  . 

We  are  now  in  a  position  to  carry  out  the  inverse  Laplace 

transform  of  r£|  (r,9,z,oj)  to  obtain   r£  (r,6,z,t).   Combining 

eqs .  (39)  and  (70)  and  assuming  eq.  (67)  for  J  (oj)  ,  we  find, 

after  deforming  the  contour  L    to  a  contour  L'    which  encircles 

the  singularities  of  the  Laplace  integrand  and  extends  into  the 

upper  half  co-plane,  the  following  expression, 

r£r(r,t)  =  J0R^(r,oj0)sinoj0t  +  JQX  (r,u>0)cosu>0t  +  n(r,t)  ,      (76) 

where  R.  and  Xf  are  respectively  the  real  and  imaginary  parts 
of  the  complex  function  Z„,  i.e., 

Z  (r,oj)  =      R_(r,w)  +  LX(r,u)    ,  (77) 


with 


Z^(r,oj)   =   rCr(r,u)/J0(u)    .  (78) 


From  eq.  (70),  we  obtain, 
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Z5(rfu)  -  uQrc  p^^-y  J(r    ,   ,  (79) 

At  r 

where  the  factor  cos (m 9  +  k  z)  has  been  dropped  for  convenience 
of  notation.   The  9  and  z  dependence  of  r^   is  then  simply 
given  in  terms  of  eq.  (76)  as, 

r5r(rf9,z,t)  =  r£r(r,t)cos(mc9+kcz)    .         (80) 

The  term  n(r,t)  is  eq.  (76)  represents  the  transient  part  of 
the  plasma  response  to  the  applied  external  force,  and  is  given 
by, 

n(r,t)  -  j  |^  exp(iut)Z  (rfu)JQ(u)  (81) 

L" 

where  the  deformed  contour  L'  encloses  only  the  singularities 

of  the  kernal  Z£(r,w)f  the  simple  poles  of  JQ  (oj)  having  already 

been  accounted  for  through  the  first  two  terms  in  the  right 

hand  side  of  eq.  (76) .   The  singularities  consist  of  the  branch 

points  and  branch  cuts  of  <Mr,u))  and  of  the  poles  of  A(w).   For 

the  profile  shown  in  Fig.  6a,  this  would  then  lead  to  the  contour 

i-'  given  in  Fig.  7  where  the  branch  cuts  have  been  analytically 

continued  off  the  real  w  axis  into  the  complex  plane.   The  poles 

indicated  in  Fig.  7  are  zeros  of  the  denominator  of  A(oj),  i.e., 

they  satisfy 

I'  (T  )K  (t  )-I  (t  )K'  (t  )  - 
m  w  m  p   m  p  m  w 


|k|qU)  [i;(xw)K;(x p>-i;<Tp >k;<tw>]  =  0  .    (82) 


The  poles  at  oj  =  toa,  w  ,  wY  are  the  magnetoacoustic  modes  which 

m   m   m 

-32- 


depend  on  the  compressibility  of  the  plasma.   For  the  incom- 
pressible plasma,  these  discrete  modes  are  absent.   The  complex 

a   6 
poles  at  cj  =  a)  ,  id  are  damped  surfaces  oscillations,  discussed 

in  detail  elsewhere  (Grossmann  and  Tataronis,  1973;  Tataronis, 

1975) ,  which  appear  explicitly  in  the  response  only  if  the  branch 

cuts  are  analytically  continued  into  the  upper  half  plane  as  in 

Fig.  7.   In  general  these  roots  of  eq.  (82)  are  complex,  and 

the  implied  decay  originates  from  phase  mixing  of  the  Alfven  and 

cusp  continuum  modes. 

Although  the  integration  of  eq.  (81)  along  L'    cannot  be 

carried  out  explicitly,  one  can  obtain  an  asymptotic  expression 

for  n(r,t)  as  t  ■*  °°.   The  result  is, 

cosoj  .  (r)  t  cosai  t 

n(r.t)%.(r,  -^♦JVD  — ^ 


+  I   S„  (r)  exp  (-u)  .  t)  cosu)   t 
a     r    si      sr 


+  I    TR(r)cosuTt   .  (83) 

.   p       m 

8 

where  w    and  cj    are  the  real  and  imaginary  parts  of  w   re- 
sr      si  s 

spectively.   The  first  two  terms  in  the  right  hand  side  of 
eq.  (83)  decay  in  time  and  arise  from  integration  about  the 
branch  points  in  Fig.  7.   The  third  term  is  the  decaying  surface 
oscillation  and  the  last  term  represents  the  steady  magneto- 
acoustic  modes. 

Our  main  interest  of  course  lies  in  regions  about  r  =  r. 

2     2         2 
where  uQ  =  w  (rQ)  or  u  (r  ) .   Here  the  singularities  of 
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Jg (W)  interact  to  yield  higher  order  singularities.  For  both  an 
Alfven  [uQ  =  ^A(r  )  ]  and  a  cusp  [ojq  =  w  (rQ)  ]  type  resonance  the 
singularity  is  of  the  form: 

£n(oj  -  oj  ) 
r?r(r,w)  ^  ^-   ,  (84) 

which  implies, 


rg    (r,t)  %  In   t  cos  u  t   ,  (85) 

asymptotically  as  t  ■+  °°.   Thus  r£  (r,t)  grows  logarithmically 
with  time  at  the  point  r  =  r  .   From  eqs.  (55)  and  (56)  it  can 
be  shown  that  ^(r/t)  ancj  £  (r,t)  at  the  point  r  =  r.  behave 
temporally  as  follows: 

CQ(r,t)  ^  5z(r,t)  ^  t  cos  coQt   .  (86) 

It  is  this  local  temporal  growth  at  rn  which  will  be  associated 
with  local  energy  absorption  at  that  point.   In  the  next  section 
we  will  use  the  above  solutions  in  an  asymptotic  calculation  of 
energy  absorption  for  t  ■>  °°. 


-34- 


5 .   Asymptotic  Form  of  Energy  Absorption 

In  this  section  we  will  consider  the  time  behavior  of  the 
plasma  energy  W   in  the  limit  as  t  ■*■  °°.   We  will  average  dW  /dt 
over  a  period  of  the  given  frequency  of  the  external  source.   If 
this  average  is  greater  than  zero,  then  it  can  be  concluded  that 
the  plasma  absorbs  energy.   Alternatively,  one  can  obtain  the 
average  rate  of  change  of  the  plasma  energy  by  calculating  the 
average  rate  at  which  energy  is  delivered  by  the  external  source. 
Because  of  eq.  (27)  which  expresses  conservation  of  energy,  the 
average  loss  or  gain  of  energy  by  the  source  should  equal  re- 
spectively the  average  gain  or  loss  of  energy  by  the  plasma. 
This  equality  will  be  established  in  Section  6. 

First  we  specialize  the  general  form  of  W  ,  the  total  perturbed 
plasma  energy,  given  by  eq.  (28),  for  the  specific  case  of  cylin- 
drical symmetry.   Assuming  periodic  boundary  conditions  in  both 
6  and  z,  a  lengthy  calculation  yields, 


dW 
dt 


d_ 
dt 


+  2 


dx 


I  V^e^z*  +  217-  {(B-V^)2+(B.VCe)2+(B.V^)2^ 


f 


rp+  5_ 


(v-c)2- 


'b2^ 


2V 


2r2   1 


2B' 

B-V(B-U+  — |  r^J 


2BA 

u„r  ^r  -    6 


(87) 


where  (  ) '  denotes  differentiation  with  respect  to  r.   The  volume 
element  dT  5  r  dr  d9  dz  with  0  <  r  <  r   ,  0  <  9  <  2tt  and 
0    z    2tt/k  .   In  eq.  (87)  we  denote  by  6KE  the  first  three  terms, 
the  perturbed  kinetic  energy  of  the  plasma.   The  remaining  terms 
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we  will  denote  by  <5w,  the  perturbed  potential  energy.   In  the 
following  we  shall  make  use  of  the  identity: 


2ir/k 


d6 


dz 


cos  (m  8+k  z)  _ 

_  =  7TA   ,   A  =  t- — 

sin  (m  9+k  z)  c 

c   c 


(88) 


The  integrands  of  <5KE  and  6W  are  functions  of  r,  m   and  k  .   We 

c       c 

have  seen  that  at  the  point  r  =  r  where  either  A  or  C  vanishes, 
singularities  of  the  Laplace  transform  of  the  fluctuating  vari- 
ables exist  with  the  consequence  that  these  variables  grow  with 
time  at  rn  upon  application  of  a  harmonic  force.  At  every  other 
point,  the  perturbation  fields  saturate  in  amplitude  as  t  ■*■<*> . 
Therefore,  for  t  large,  the  major  contribution  to  the  kinetic 
and  potential  energies  arises  from  some  narrow  region  of  width 
Ar  about  the  point  rfi.   We  express  this  as: 


6KE  = 


6W  =  ^— 


rQ+Ar 


rQ-Ar 


ttX 


rQ2 


rQ+Ar 


rQ-Ar 


r  dr  pv   +  6  (6KE) 


(89) 


r  dr 


=—  +   Q-(JxC)  +  V-US  P'  +  rpv-o 
M0 


+6  (6W)  (90) 


Using  here  again,  the  general  form  of  eq.  (28),  section  2. 
The  terms  6 (6KE)  and  6 (6W)  result  from  the  remaining  parts  of 
the  spatial  integration.   As  t  +  »,  d6(6KE)/dt  and  d5(6W)/dt  •*■   0. 
Therefore  in  computing  the  time  rate  of  change  of  the  energy 
6(6KE)  and  6(6W)  can  be  neglected.   Considering  the  solutions 
found  in  the  last  section  we  notice  that  the  fastest  growing 
fields  and  variables  at  a  singular  point  r_  are  v  ,  v   ,  £„  ,  £   , 
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QQ  /  Q    for  both  the  Alfv£n  and  the  cusp  mode   with  the  excep- 
ts    z 

tion  that   V«£   must  be  included  for  the  cusp  singularity.  In 

the  following  we  treat  the  Alfv6n  and   cusp  absorption  separately. 


Alfv£n  Absorption 

2     2 
At  the   point   rQ   where  ajQ  =  w  (rQ),  the  function  A(rQ,uj) 

vanishes,  and  it  can  be  shown  that  Qa    ,  Q   ,  £Q  /  5   /  vD  and 

0       Z      0       Z      o 

v   are  related  accordingly, 


QQ  =  F  5 


9  ' 


Q   =  F  £   , 
z     ^z 

2     2   2     2 

v2  -  v2  +  v2  =.2(rQ)(C2  +  #  , 


(91) 
(92) 
(93) 
(94) 


where  the  definition   y  =  9§/3t  has  been  used  in  eq .  (94). 
Consequently,  the  form  of  6KE   and  <5W  reduce  to  the  following 
expressions , 


irXr, 


6KE 


tt  Ar, 


6W  = 


ro+Ar 


rQ-Ar 
rQ  +  Ar 


r0-Ar 


dr  pv 


(95) 


dr 


(96) 


where  for  clarity  we  have  neglected  6(6KE)  and  <5  («5W)  .   Substituting 
ecl-     (93)  into  eq.  (96)  we  obtain  an  expression  for  the  lowest 
order  total  perturbed  plasma  energy  about  the  singular  point, 
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TiAr, 


SKE    +    Sw 


P(r  ) 


vAr 

f 


dr 


rQ-Ar 


2  2  2  2       2   ^ 

.Ve  +  vz  +  VV(WJ 


(97) 


2         2  • 

where  wA(r0)  =  F  (rQ) /p(rQ) pQ.   In  order  to  proceed  it  is  neces- 
sary to  obtain  the  variables  vfl(r,aj),  v  (r,u>),  £fl(r,a))  and 


?  (r,ii))  about  the  point  rQ, 


This  can  be  readily  accomplished 


since  the  form  of  r£  (r,ui)  at  this  point  is  known  and  is  given 
by  eq.  (84) •   Using  eqs.  (55)  and  (56)  we  can  define  £Q   and  Cz 
in  terms  of  r£   and  (r£  ) ' ;  the  inverse  Laplace  transform  may 
be  accomplished  with  the  result  in  the  limit  t  -*-  °°, 

,  o  o   Rl+Rl         ■>  sin2[l  aJA(rn)  (r_ro)t] 

d(r,t)    +   ej(r,t)    =   Jl   -^  cos2(^t+*J 


4u>, 


"0 


^^(ro)(r-ro)]2 


(98) 


and 

.  -  ,   R2+R2  j  sin2[|  w'(r  )(r-rfl)t] 

v2(r,t)  +  v2(r,t)  =  J2  -^  sin2(a)nt+4)fl) 


4a), 


[|^(r0)(r-r0)]2 


(99) 


where  RQ  and  R  are  constants  which  are  specified  by  applying 

D  Z 

boundary  conditions  at  the  plasma  vacuum  interface  in  the  same 
manner  in  which  C(w)  was  determined  for  the  function  r£r(r,ui) 
through  eq.  (60).   Using  eqs.  (98)   and  (99)   we  can  express  the 
total  energy  as, 


6KE    +    6w   = 


•nXr 


2      2    r0+Ar 
0       ,      .    2,       ,     T2    R6+Rz 
2-  p(r0)uJA(r0)    J« 


4oj 


dr 


0 


r0-Ar 


sin2[|  ^(rQ)  (r-rQ)t] 
[|^(r0)(r-r0)]2 


(100) 
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irAr, 


r2  Re+Rz 
re  2~~ 


2TTt 


(101) 


Applying  the  boundary  conditions  on  £fl  and  £   through  the 
corresponding  boundary  conditions  for   r£   yields  the  result: 


*l*l  - 


where 


4      2    2 
F*    BE 


LHo 


roD 

0        J  r. 


rX<vi2 


2    2 

U0U0      _ 

**£">!  W  -k<<yi2 

x    [<"l(r0))'I2 
l<MO  I2 


(102; 


D(rQ)    - 


F2(rQ) 


E2(rQ) 


(103) 


and  in  the  above,  for  the  purpose  of  clarity,  the  wall  radius 
r   has  been  removed  to  infinity.   Substitution  of  eq.  (102) 
into  eq.  (101)  yields  the  final  result  for  the  temporal 
behavior  of  the  total  perturbed  plasma  energy  in  the  singular 
layer, 


6KE  +  5W 


^    2T2 
*XrcJ9 
4r„ 


(nunt) 


B2(rQ) 


ETTr^7  F2(rp) 


p(rQ) 


m   c 


IVTp}  "  ^  Km(V' 


|(^(r0))  | 


(104) 


We  now  differentiate  (104)  with  respect  to  time  to  obtain  the 
power  supplied  to  the  plasma  and  average  this  expressicn  over 
one  period,   T.  ,  of  the  external  source,  T-  *  2~/ui      ,    to  obtain 
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^c       T2         B2(r0}    4P{T0] 


[Km(Tc)]2                     l<ui<ro,,,|2 
x  HL-E a_u (105) 

Equation  (105)  is  a  positive  definite  quantity  representing  the 
energy  absorbed  by  the  plasma  in  a  narrow  region  about  the  point 

rfl  where  the  external  frequencey  of  the  source   oj-  equals  the 

2 

local  Alfv£n  frequency  of  the  plasma  (joA(rQ). 


Cusp  Absorption 

2     2 

We  now  consider  the  case  when  coQ  =  uc(rQ)  .      The  procedure 

here  is  identical  to  that  for  the  Alfven  mode.   As  mentioned 
earlier,  the  fastest  growing  fields  at  the  point  r_  where  C 
vanishes  are   QQ  ,  Qz  ,  vQ  ,  v^    ,  EQ  ,  Cz  and  V«E.   The 
divergence  term,  V*£  ,  reflects  the  fact  that  the  cusp  mode 
only  occurs  when  compressibility  is  included  in  the  equations 
of  motion.   Again  using  the  relationships  between  £e  ,   E 
and  rE   it  can  be  shown  that  where  C(rQ)  vanishes, 

QQ  =  FCe  -  BqV-C  ,  (106) 

Q   S  FE   -  B  V-E  .  (107) 

The  divergence  V»E   can  similarly  be  treated  and  at  r.  can  be 

shown  to  be 

A(rn}  h 

v.£  = °— _  |_  rE    .  (108) 

-?         ,   .  2  dr   ^r 
r0p(r0)o) 
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The  form  of  the  total  perturbed  energy  for  the  cusp  mode  is 
given  in  a  similar  manner  to  that  of  the  Alfven  mode,  i.e., 
we  consider  again  only  the  contribution  arising  from  some 
narrow  region  of  width   2Ar   about  the  resonance  or  singular 
point   rQ   and  neglect  the  contributions  from  the  remainder 
of  the  integration.   We  have  then  for  the  cusp  mode: 


itXr, 


rQ+Ar 


6KE  = 


dr  p(r)  [v?(r)  +  v* (r) ]  , 
o        z 


(109) 


rQ-Ar 


frXr, 


rQ+Ar 


<5W  = 


dr 


rQ!?(r)    Q2(r) 
— +  -^ +  rp(r)  (v^) 


rQ-Ar 


.  (110) 


2         2 
We  now  form  the  products   QQ   and   Q   , 


Q2  =  F2Cg  +  Bq(V-C)2  -  2BQF5Q  V-C  , 


Q2  =  F2^2  +  B2(V-C)2  -  2B  F^  V-E,    , 
z      z     z    -       z   z 


(111) 

(112) 


such  that, 


Q2  +  Q2  =  F2(Se+S2)  +  B2(V-C)2  -  2FB-§  V-| 


(113) 


Using  once  more  the  relationship  between   £Q,  t,      and   r^    we 
find  at  the  point   rfl   where   C(rQ)  =  °  ' 


LSL  {-TPA  E  +  pJ    E  %    ±-   r£ 
*6    rD       r         vQ     dr   sr 

n 

•-,  =   7H  (rpAk  +  pw2e  r1}  T~  r 

z    rD  uQ      dr    r 


(114) 
(115) 
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where 


D(rQ)  =  P2o>4 


(116) 


Expanding  the  numerators  of  eqs .  (114)  and  (115)  leads  to  the 
simpler  expressions, 

2  2    ^6F 

-    TPAm/r   +    pu  B  E/y_   =   -    pu>     —— 

Z  U  M/i 


2  2    BzF 

TPAk    +  PUJ    BfiE/yrt    =  PO)      -7— 

6    0  P0 


(117) 
(118) 


Hence,  the  form  of   CQ   and   C    can  be  conveniently  written  as, 


2  dr  r*-r  ' 


yQrpw 


2  dr  r^r  ' 


y0rpu 


from  which  we  obtain, 


?5  „  2   dr   ^r 

yQrpw 

2   2 

An  expression  for   5Q+?    can  now  be  obtained  as, 


(119) 


(120; 


(121) 


<i*<2 


B2F2 


2    2    2    4 

r   yQp   a) 


dTr5r 


(122) 


from  which  the  following  is  obtained 


2       2       2 

(b>b:)f*a   rd 


<§'£>F    V>§    ~~         2    2    2    4       [d?r5r 


2       2       2 
F2(^H2} 


r   yQp    a) 

B2F4 

2    2    2    4 

r  yQp   co 


3F«r 


(123) 


(124) 
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It  is  convenient  to  consider  the  terms  6KE  and  6W  separately; 
we  begin  with  6W.   So  far  we  have  obtained  the  following  form 
for  6w, 


irXr, 


6W  = 


rQ  +  Ar 


dr 


rQ-Ar 


B2p4      yP+B2/M0)A2    2B2F2A2 


2    3    2    4 

r  uqp  u 


2    2    4 

r    p   a) 


2    2    4   5 

UQP    oj  r^ 


[fert, 


(125) 


which  is  equivalent  to, 


5W  = 


2r„ 


2  4 

p  co 


rQ+Ar 

f    / 
dr 

rQ-Ar 


I?  r5r 


(126) 


A  computation  of  the  inverse  Laplace  transform  of  dr£  /dr  yields 
the  result, 


fe*«r«« 


2w 


y  |R(u>0)  I  2  cos  (u>0t-H|>) 


sin[?^(r0)(r-r0)t] 
tu»Mr0)(r-r0)] 


'  (127) 


where  R(u>n)  is  a  constant  which  will  be  specified  later.   The 
integral  in  eq.  (126)  can  thus  be  carried  out  to  give  the  result 


rQ+Ar 


dr 


I?*, 


r0-4r 


2  2 

Jfl  ~>     COS   (u)_t  +  (Jl)TTt 

-|  |RU0H2  2^ 


(128) 


An  evaluation  of  boundary  conditions  at  the  plasma-vacuum  inter- 
face allows  R(wn)  to  be  specified.   The  result  is, 


2  2  2 


IK  •  <  TT„>  ] 

m   c 


[C-J(r0))'l2 


pr<v  iKm(v  -  kqK;(v|2    i va*s i 2 


(129) 


-43- 


where  again  for  clarity  the  wall  position,  r  =  r   ,  has  been 
taken  to  be  at  infinity.  Expanding  the  coefficient  before  the 
integral  in  eq.  (126)  yields  a  simpler  expression  which  when 
combined  with  the  result  of  the  integration  yields  the  final 
form   for   6W 


TtXr, 


6W   = 


2    2 
B   F 


Jq|R(w0)  I2    irt   cos2(w0t+<J>) 


2  2  4  2, 
^0P  u0r  r=r, 


to, 


|(^(r0))'| 


(130) 


The  calculation  of  the  kinetic  energy  proceeds  in  a  similar 
manner.   We  first  express  vQ  and  v  as  functions  of  t  through 
the  relations  given  by  eqs .  (55)  and  (56).   The  result  is 


ve(t)  = 
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|R(«V 
2~"  Je   2ZI7" 

rQU0PU)0  0 


sin[^  ^(rg)  (r-rQ) t] 


2  sin  (wnt+4))  — 
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(131) 


and 


v2(t)  =  - 


B  F 
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|R("0) 

~ 2  Je   217T" 
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2  sin  (u)0t+<}>) 


sinfj  uc(r0)  (r-ro)t^ 


[<^(r0)(r-ron 


(132) 


where  we  have  retained  only  the  parts  which  lead  to  growth 

in  5KE  at  r  =  rQ.   Collecting  results,  the   expression  for  6KE 

at  this  point  can  be  written  as, 


<5KE   =  —T—  p(rQ)       nu 


|R("0)  I 


y0p   u0r      r=r0  4U)0 


2 
sin    ((jj_t+4>) 


rotAr 
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[^'(r0)(r-r0)]2    ' 


(133) 
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which  integrates  to  yield, 


trXr, 


6KE  = 


p(rQ) 


B2F2 


2  2  4  2. 

0    0    r=r, 


2  |R(u0) 


x  sin   (a)  t  +  <J>) 


IT  t 


|(<V(rn))'! 


(134) 


If  we  now  add  6KE  to  SW  we  obtain  the  final  result  for  the 
total  perturbed  plasma  energy, 


TT2Xr2        2  B2(r  ) 
6KE  +  6W  =  — ^S.  U(Jty0  J^  — ^ 


r0E  V 


i*i,  1 
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JcVio 


[K  (t  )] 
m  c 


'W  "  k^(V|2   !♦(*«>  I'  ' 


(135) 


2      2 
where  here  u   -  (F  /pun)  _    •  Taking  the  derivative  with  respect 
A         u  r-rQ 

to  time  of  eq .  (135)  gives  the  following  result  for  the  power 
supplied  to  the  plasma, 


m   c 


IVV  -  ^VV1    l*(rpH 


(136) 


Equation  (136)  is  a  positive  quantity,  and  we  interpret  this  as 
a  demonstration   that  energy  emanating  from  an  external  source 
may  be  absorbed  by  the  plasma  at  a  point  r0  where  the  external 
exciting  frequency   uQ   equals  the  local  cusp  frequency   u  (rQ) . 
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We  have  now  demonstrated  that  it  is  possible  through  wave 
resonance  phenomena  at  the  local  Alfven  and  cusp  frequencies 
in  an  inhomogeneous  plasma  column  to  couple  energy  from  an 
external  source  to  the  plasma  in  the  resonance  layers. 

Even  though  the  expressions  for  the  absorbed  power  are  convin- 
cing,  a  more  useful  quantity  for  the  purposes  of  numerical 
computations  is  the  effective  coil  impedance  which  is  felt  by 
the  source  due  to  the  presence  of  the  plasma.  In  the  next 
section  we  calculate  this  impedance  and  use  it  later  for 
numerical  investigations  of  energy  absorption  in  several  real- 
istic plasma  configurations. 
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6 .   Coil  Impedance 

A  very  useful  and  measurable  quantity  which  indicates  the 

amount  of  absorbed  and  circulating  energy  in  a  system  which  is 

excited  by  an  external  coil  is  the  coil  impedance.   Usually  the 

impedance  is  complex,  the  real  part  indicating  absorption  (or  in 

classical  systems  dissipative  damping) ,  and  the  imaginary  part 

representing  the  fraction  of  the  total  energy  circulating  in  the 

system.   We  will  see  that  the  average  charge  of  the  plasma  energy 

<dW  /dt>  as  computed  with  the  resistive  part  of  the  coil 
P 

impedance  is  identical  with  the  results  of  Section  5  where 
<dW  /dt>  is  computed  directly.   The  current  in  the  idealized  coil 
which  we  have  chosen  flows  on  the  coil  surface.   We  begin  the 
calculation  of  the  impedance  by  noting  that  the  time-averaged 
power  absorbed  by  the  plasma  is  obtained  from  eq.  (32)  and  is 
given  by  the  real  part  of, 

P  =  -  ^  ttL  r   6E-J*  ,  (137) 

2    z  c   -  -s 

where  P  is  the  absorbed  power,  6E,  as  defined  before,  is  the 

* 
perturbation  electric  field  at  the  coil  and  J   is  the  complex 

conjugate  of  the  total  coil  current  which  has  components  in  the 

9  and  z  directions,  as  given  by  eq.  (5).   Use  of  Maxwell's  equation 

plus  eq.  (9)  allows  5E'J   to  be  written  as 


• 


u 


5E-J    =   -  =•  &B     .  J     ,  (138) 


•  s        k    vr 
and  hence  the  total  time-averaged  power  as, 


P  -  4  nL  r   £  6B    J  (139) 

2    z  c  k    vr   so 


6b    is  derived  from  the  magnetic  potential  <j>  where  6B  =  V<j>  and 
$   is  given  by  eq.  (A8)  in  the  Appendix.   Computing  6B   as  9<j>/9r 
and  substituting  into  eq.  (139)  yields  the  following  equation 
for  the  average  power  P, 

7rrc  rc  X2U)  ,    ,2 
P-  -  lu>yQ  -^-E—-  |jse|  (140) 

where  x?^)  anc^  A(ui)  are  given  by  eqs.  (64)  and  (70b)  respectively. 
In  eq.  (140)  k  and  m  are  respectively  equal  to  m   and  (-k  ) . 
From  this  expression  we  can  define  an  impedance  Z(uj)  given  by: 

P  -  j  Z(w) |Jse|2  ,  (141) 

where 

r   x2(u) 

Z(uj)  =  -  icoMrt  irr   —  .  ,    ,  (142) 

0    c  r    A(u>) 
P 

Z(cj)  will  have  both  real  and  imaginary  parts  if  q  is  complex. 
The  value  of  q  is  obtained  from  knowledge  of  the  radial  component 
£   at  the  plasma  boundary.   Calculation  of  5   and  the  subsequent 
value  of  q  is  discussed  later.   It  is  clear  however  that  in  the 
absence  of  singularities  in  the  Alfven  or  cusp  continuum,  £   is  a 
real  quantity  resulting  in  a  real  value  for  q,  and  hence  the 
impedance  Z(co)  will  be  purely  imaginary.   This  would  imply  no 
absorption  and  the  total  power  in  the  system  would  be  circulatory. 
We  point  out  that  the  real  part  of  eq.  (141)  is  identically  equal 
to  eq.  (105)  in  the  Alfven  continuum. and  to  eq.  (136)  in  the  cusp 
continuum.   Therefore,  the  absorbed  energy  implied  by  the  real 
part  of  Z(w)  is  transferred  to  the  local  resonances  of  ideal  MHD. 
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7 .   The  Plasma  Profiles 

We  present  here  the  actual  plasma  profiles  which  have  been  used 
to  calculate  the  absorption.   As  shown  by  Tataronis  (1975)   the 
quantitative  absorption  rate  depends  sensitively  on  the 
functional  shape  of  the  continuum  of  Alfv6n  frequencies  inside  the 
plasma  column.   This  means  that  although  positive  absorption  can 
be  found  for  any  plasma,  significant  differences  may  be  noticed 
between  several  different  plasma  configurations.  By  using  real- 
istic profiles  we  are  assured  that  any  peculiar  characteristics 
of  real  experimental  plasmas  are  taken  into  account.  For  our 
calculations  we  have  chosen  both  Tokamak  and  Scyllac  profiles. 
The  profiles  have  been  assembled  through  available  experimental 
data  and  in  some  cases  best  guesses  averaged  over  several 
experiments.   It  should  be  pointed  out  here  that  the  profiles 
have  been  adjusted  to  exclude  effects  of  toroidicity  since  the 
present  calculations  are  carried  out  for  straight  cylindrical 
geometry . 

Details  of  the   experimental  measurements  are  not  discussed 
here.   The  first  case  presented  here  is  the  Princeton  ST  Tokamak. 
This  plasma  is  basically  a  low-8  plasma  with  the  equilibrium 
being  given  mainly  by  the  poloidal  magnetic  field.   The  functional 
form  for  the  various  profiles  are  given  below: 

B  (r)  -  BQ  -  const.  ,  (143) 

p(r)  -  P0(l  -(r/Rp)4l ,  (144) 
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|    (1    "     <"^>4) 
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1/2 


,    r    >  2  ,  -1 ,    r    . 2 

(— -)    +  sin       (-7—) 

P  P 


(146) 


In  the  above  pn   is  the  density  on  the  axis,  I   is  the  total 
longitudinal  (toroidal)  current  in  amperes,   r  the  plasma  radius. 


Equation 


evaluated  at  the  plasma  radius  yields , 


«  1   1    y°Iz 
B9(rp)  =  2ir" 


(147) 


The  equilibrium  relationship,  eq.  (4) ,  gives  the  radial  pressure 


variation  as : 


2„2    n 


P(r)    = 


u  b 


'0 


l-f 


L+  f  fe(x')2  dx* 
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-I 


;   fe(x')2   dx' 


x' 


(148) 


where    fQ   =  BQ(r)/(yB.)    and   the   parameter   y   is   the   ratio  of   the 
WW  u 

poloidal   to   the   longitudinal     magnetic   field  at   the  plasma 
radius,    i.e. , 


V    = 


_Lofz 

27TrPBo 


(149) 


In  eq.  (148),  x  =  r/r  .    For  these  profiles  the  limit  of  P/p 

XT 

as  r  ■+  r   yields  zero.   At  the  plasma  radius  r  =  r   ,  vacuum 

p  *  f  p 

profiles  for  Bfl  are  joined,  and  P  and  p  are  identically  zero. 
The  following  table  contains  values  of  the  constants  in  the  above 
formulas  obtained  from  the  experimental  data  of  the  ST  Tokamak 
(Hosea,  1974) : 
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Parameters  of  the  ST  Tokamak 

BQ   =   12  -  45  k  Gauss 

pQ   =   2  x  io13  /cm3 

I    =   10  -  100  k  amperes 

r   =   12  cm 
P 

The  distribution  of  Alfv^n  frequencies  in  the  plasma  column  is 
calculated  from  the  following  equation  after  the  parameters  "m" 
and  "k"  have  been  chosen: 

2 


1    T  m 
>(r)  L  x 


"a'^TT  |-^fe(x)  +k 


(150) 


where  u.  is  the  frequency  normalized  according  to 


id  =  1-  ,     v   =  ,   and    k  =  kr 

P 


Similarly,  the  cusp  frequency  is  given  by  the  expression, 
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(TP+l+U   fg) 


T-  uffi  +  kl   '  (151) 


2 
where   P  is  the  normalized  pressure,  (P/BQ/uQ) .      Figure   8 

presents  the  radial  dependence  of  the  profiles  of  B   ,  f.  ,  P 

and  p  from  the  axis  to  the  plasma  radius.   Figure  9   shows  the 

radial  dependence  of  oj.  and  w   .   The  value  of  y  used  for 

A        C 

Figures  8  and  9   has  been  set  equal  to  0.05.   The  values  of  m 
and  k   are  1  and  1.1   respectively.   The  value  of  k  corresponds 
to  a  60  cm  wavelength.   Near  the  edges  of  the  plasma  the  Alfv£n 
continuum  tends  to  large  frequency  values  whereas  the  cusp 
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continuum  goes  to  zero.   This  is  a  consequence  of  the  limit 

of  P/p   tending  to  zero.   From  Figure  9   one   sees  a  continuous 

range  of  frequencies  from  zero   to  infinity  with  the  exception 

of  the  strip  oi      <  u)  <  w.  It  is  readily  seen  that  by 

c      c  max  —   —  A  min  J  J 

adjusting  the  "m"  and  "k"  for  a  given  profile  this  strip  can  be 
made  as  small  as  possible. 

The  second  set  of  profiles  investigated  are  those  from  the 
ORMAK  experiment  at  Oak  Ridge  National  Laboratory.  These  profiles 
correspond  to  two  operating  regimes  of  the  experiment,  the 
so-called  "A"  and  "B"  type  of  discharges.   It  will  be  convenient 
to  present  these  profiles  in  the  following  manner.  Let  j^-(r) 
denote  a  particular  radial  profile  and  consider  ^"(r)  to  be 
given  by  the  following  function: 

m 


J-ir)    =  >(0)  (1  -  (r/rp)n)   ,  (152) 

/~> 
where  d~(0)  represents  the  function  evaluated  on  the  axis.   The 

following  table  (England,  1975)  then  gives  the  values  of  n,  m 

and  J~iO)      for  the  cases  treated  here. 

Parametres  of  the  ORMAK  Tokamak 


J-(r) 

Discharge 

Type 

n 

m 

^(0) 

P(r) 

A 

2 

2 

13    3 
4  *  10   /cm 

B 

2 

2 

3  x  I013/cm3 

Jz(r) 

A 

1.5 

1 

2 

13  3  amp/cm 

B 

4 

3 

2 

125  amp/cm 
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The  plasma  radius  is  2  3  cm  and  the  main  longitudinal  magnetic 
field  was  18  k  Gauss.   Again,  here  the   longitudinal  magnetic 
field  is  constant  and  the  poloidal  field  is  computed  from 
Ampere's  law  as  before.   The  pressure  is  given  by  eq.  (148) 
once  the  poloidal  field  is  known.   Although  we  do  not  show  them, 
curves  similar  to  those  of  Figure  9  can  be  obtained  from  the 
Ormak  profiles.   Although  similar,  the  shapes  are  different; 
this  difference  being  a  key  factor  in  the  actual  values  of 
computed  absorption. 

In  an  attempt  to  connect  the  parameters  "k"  and  u  with 
experimental   configurations,  let  us  form  the  following  dimension- 
less  number   kr  /u .   If  the  true  toroidal  plasma  had  a  circumfer- 
ence of  length  L  =  27:1^,  where  FL,  is  the  major  torus  radius,  then, 

kr     _     B  (r  )  B-.r 

where   n   is  the  number  of  wave  lengths  along  the  cylinder  of 
length   L  and   q(r  )   is  the  well  known  safety  factor  evaluated 
at  the  limiter.   Typically,  for  the  range  of  parameters  examined 
in  this  work,  the  value  of  q  at  the  plasma   radius  (limiter  in 
the  experiment)  was  3  and  n  corresponds  to  the  order  of  about 
10  wave  lengths. 

We  turn  next  to  the  high-8  Scyllac  profiles.  Although  the 
Scyllac  experiment   is  basically  a  helically  toroidal,  large 
aspect  ratio  (R  /r   >>  1)   plasma  configuration,  the  main  magnetic 
field  component  is  the  6-pinch  field.   All  other  fields  are  of 
a  much  smaller  order.  For  this  reason  we  will  approximate  t I 
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Scyllac  plasma  by  a  0-pinch.   We  mention  only  briefly   that  the 
existing  Scyllac  feedback  stabilization  system  could  be  effec- 
tively used  to  investigate  the  Alfv£n  wave  absorption  (assuming 
stable  plasmas  with  sufficient  life  times  can  be  obtained) . 
Here  there  is  only  one  field  component,  the  axial  field,  and 
we  use  as  an  approximation  the  following  function: 

Bz(r)  =  -  (BQ-Bi)  [(r/rp)4  -  2(r/rp)2]  +  Bj_  ,  (154) 

for  0  <_  r  <_  r   ,  and, 

B   =  Brt     for     r  >  r   .  (155) 

z     0  p 

B_  is  the  constant  B    field  in  the  vacuum  and  B.  is  the  value 
0  z  i 

of  B  at  r  =  0 .   The  standard  definition  of  B  is  thus, 
z 

f  21 1/2 

■  B  -  [l  -  (Bi/B0)2j 

From  eq.  (4)  the  pressure  is  found  as, 

»     B°  f,     B*  1 

P  "  W0    I1  "  ^  J  <"« 

For  the  density,  a  Gaussian  profile  adjusted  to  zero  density  at 
r  =  r    and  further  adjusted  for  an  appropriate  half  width 
consistent  with  the  above  equilibrium  is  chosen.  Typically  the 
magnetic  field  strength  B_  was  taken  as  40  k  Gauss,  the  density 
on  axis  to  be  3  *  10   /cm  ,  the  plasma  radius  2.5  cm   and  the  wall 
radius  5  -  10  cm.   The  value  of  B  can  be  varied  from  zero  to  unity. 
Figure  10   shows  the  radial  variation  of  the  equilibrium  quantities. 
The  Alfve*n   and  cusp  continuum  frequencies  are  shown  in  Figure  11 
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for  the  choice   m  =  1  and  k  =  0.4.   The  value  of  B  as  shown  in 
the  curves  was  taken  to  be   8  =  0.5   on  the  plasma  axis. 

We  mention  that  the  profile  given  in  eq.  (155)  for  the  B 
field  can  be  used  to  construct  reversed  field  conditions  such 
as  those  produced  in  a  reverse  bias  6-pinch  experiment.  For 
this  situation  the  8  definition  must  be  altered  to, 

B  =   max   ,  (157) 

Bo/:% 

where   P    is  the  value  of  the   pressure  at  the  point  where 

max  c 

the  B   field  reverses.   For  a   8  =  0.5   we  show  in  Figure  12 
the  Alfve'n  and  cusp  continuum  again  for  m  =  1  and  k  =  0.4. 
Notice  that  it  is  now  possible  to  find  multiple  points    in  the 
plasma  where  the  applied  frequency  equals  a  local  Alfv6n  or 
cusp  frequency.   It  has  been  conjectured  that  such  profiles 
lead  to  enhanced  absorption  (Tataronis,  1974).   The  plasma 
produced  in  the  Garching  poloidal  belt-pinch  experiment   is 
an  example  of  such  a  configuration. 

We  next  discuss  certain  numerical  techniques  used  in  applying 
these  profiles  to  the  energy  absorption  formulation. 
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8 .   Numerical  Techniques 

In  this  section  the  pertinent  numerical  techniques  used  in 
obtaining  the  final  value  of  the  coil  impedance  are  discussed 
briefly.   As  shown  in  Section  6,   the  impedance  calculation 
requires  the  value  of  the  quantity  (r£  )  '/(rE,    )   evaluated  at 
the  plasma- vacuum  interface.   This  quantity  is  obtained  as 
follows.   The  governing  equations  (40)  and  (41)  for  rC   and  p. 
respectively  are  expanded  about  the  origin  to  give  starting 
conditions  for  the  numerical  integration  of  these  two  equations. 
A  convenient  normalization  specifies  the  initial  data.   From 
the  origin  up  to  the  position  in  radius  where  the  first  Alfven 
or  cusp  singularity  occurs,  the  solutions  for  r£   and  p,  are 
purely  real.   The  treatment  of  the   solutions  about  singulari- 
ties are  handled  by  defining  a  small  region  of  width  2e  about 
the  singular  point.   The  differential  equations  defining  r£ 
and  p*   are  integrated  numerically  up  to  the  point  rQ-e  ,  where 
rQ  denotes  the  singular  point.   These  solutions  are   then 
continued  anaytically  to    the  other  side  of  the  singularity 
at  r  =  rn+e  by  means  of  a  power  series  expansion  of  the  solution 
as  given  by  eqs .  (52)  and  (53).   Specifically  we  invoke  for  oj 
real  the  following  analytic  continuation, 

a,  (u)  [a-  (w)g(r,w)+g(r,u>)  In  (rn-r)  +  (r-rn)  h  (r  ,u>)  ]  ,  r<r_,  (158) 


sr 


a1(u)  [  (a2+iir)g(r,cj)  +g(r/u)  In  (r-rQ)  +  (r-rQ)  h  (r  ,u>)  ]  , 


>rQ,  (159) 


where  the  constants  a,  (id)  and  a2(w)  are  determined  by  assuming 
continuity  of   r£   and  dr£  /dr   at  r  =  rQ-e,  and  the  sign  of  i 
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in  eq.  (159)  is  uniquely  determined  by  invoking  causality  in 

the  form  of  the  limiting  process    Im(uj)  -*-  0  from  the  lower 

half  of  the  complex  ui  plane.   At   r  =  rn+e,  the  values  of  r£ 

and  p,  serve  as  initial   data  in  order  to  continue  the  numerical 

integration  to  r  =  r  .   Up  to  the  position  of  the  first 

singularity,  the  solutions  are  purely  real.   Continuation  past 

the  singularity  with  eq.  (159)  results  in  an  imaginary  part 

from  the  logarithm,  and  therefore  the  solutions  will  be  complex 

at  r  =  r  .    At  this  point  the  value  of  q,  eqs .  (60)  and  (65) ,  can  be 

calculated  and  hence  the  absorption  rate.   The  dependence  of 

the  value  of  q  on  the  width  of  the  region  2e  was  investigated  by 

varying  both  the  numerical  size  of  e    and  the  number   of  terms  in 

the  series  expansions  for  r£   and  p.  about  the  singular  points. 

Both  e  and  the  number  of  terms  were  varied  until  convergence 

on  q  was  obtained. 

One  special  case  must  be  mentioned.   It  is  possible  that  for 
certain  profiles  at  a  given  frequency,  the   position  of  a  singular 
point  may  fall  very  close  to  a  zero  of  the  fast  and  slow  magneto- 
acoustic  branch.   When  two  such  points  identically  coincide  the 
singularity  disappears  (the  strength  is  reduced  from  In  (x-x_) 
to  (x-xn)  In  (x-x  ) ) .   A  problem  in  the  numerical  treatment  of 
such  a  small  separation  of  singularity  and  zero  arises  due  to 
the  occurrence  of  large  gradients  in  the  region  between  the 
points.   This  can  be  eliminated  by  using  an  inner  and  outer 
expansion  about  the  singularity  and  zero  respectively  and  analyti- 
cally continuing  past  the  mutual  region  where  both  expansions  are 
valid.   For  the  cases  treated  by  us,  such  problems  only  occurred 
for  frequencies  in  the  cusp  continuum. 
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All  numerical  work  was  carried  out  on  the  New  York  University, 
Courant  Institute  CDC  6600  Computer.   Typically  a  computation  of 
the  impedance  for  one   frequency,  one  wavelength   and  "m"  value 
requires  about  one  second . 
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9 .   Numerical  Results 

In  this  section  we  present  the  results  of  numerical 
computation  of  the  effective  coil  impedance  for  the  profiles 
specified  in  Section  7.  The  pertinent   numerical  details 
concerning  the  integration  of  the  equations  of  motion  (40) 
and  (41)  have  already  been  given  in  Section  8.   We  begin  with 
a  discussion  of  the  ST  Tokamak  plasma. 

The  available  frequency  range  for  both  the  Alfven  and 
cusp  bands  can  be  computed  from  eqs.  (150)  and  (151). 

From  Figure  9   an  indication  can  be  obtained  as  to  the  position 
where  the  iraximum  absorption  for  the  Alfven  waves  occur  .  This 

is  true  since  the  power  absorbed  is  directly  proportional  to 

2 

|  (u)A(r)  )  '  |   at  the  resonance  point  as  shown  by  eq.  (104)  . 

From  Figure  9   it  can  be  seen  that  strong  absorption  should 
only  occur  at  radial  positions  greater  than  at  least  60%  of  the 
plasma  radius.   This  will  be  shown  to  actually  be  the  case  as 
borne  out  by  the  numerical  results.   Figure  13   presents  the 
results  of  the  numerical  computations   of  the  real  part  of 
the  coil  impedance,  Z   ,   versus  the  frequency  f.   Both  Z   and  f 
are  given  in  dimensional  form,  ohms  and  cycles  per  second  respec- 
tively.  The  first  noticeable  result  is  the  fact  that  the  cusp 
impedance  values  are  very  small,  typically  less  than  about  0.02 
ohms  over  the  cusp  frequency  range.   For  this  reason  they  are 
not  shown  plotted  in  the  figure.    The  Alfve'n   resistance  however 
is  plotted.   The  values  for  the,  admittedly  very  ideal  case,  are 
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impressively  large.   Here  the  ratio  of  coil  radius  to  plasma 
radius  is  1.5  and  the  wall  has  been  taken  at  infinity.  We  will 
discuss  in  detail  later  the  important  effect  of  the  conducting 
wall.   Shown  in  the  figure   are  the  results  for  both  an  m  =  1 
and  m  =  2  coil.   The  effect  of   changing  m  for  the  ST  plasma 
is  seen  to  be  rather  slight,  the  main  effect  being  a  shift  in 
the  position  of  maximum  absorption.  For  the  ST  data  the  absorp- 
tion curves  rise  sharply  to  a  maximum  and  then  fall  gradually; 
the  range  of   frequency  for  strong  absorption  is  rather  broad. 
No  attempt  was  made  to  try  to  maximize  the  impedance  by  varying 
the  coil  wave  length;  the  value  was  chosen  to  be  60  cm  and 
corresponds  to  about  10  wave  lengths  around  the  actual  toroidal 
plasma.   It  is  also  noticed  that  the  absorption  falls  to  zero 
at  isolated  values  of  frequency.   Analysis  of  this  behavior  shows 
that  the  vanishing  of  the  resistance  at  isolated  frequencies 
arises  if  the  solution  of  the  MHD  equations  are  simultaneously 
analytic  at  the  Alfven  resonance  and  consistent  with  the  boundary 
condition  at  r  =  0.   One  then  concludes  that  the  logarithm 
appearing  in  eq.  (5  3)  is  absent,  implying  that  £   and  q  are  real, 
and  hence  that  Z  is  imaginary.   An  inspection   of   5r(r)  for  <j  in 
the  neighborhood  where  the  phenomenon  occurs  shows  that  £r 
oscillates  in  r,  with  one-half  period  of  spatial  oscillation 
entering  the  domain  between  r  =  0  and  the  singular  point  at  each 
crossing  of  a  critical  to.   These  oscillations  probably   corres- 
pond to  the  excitation  of  the  fast  MHD  wave  which  propagates  in 
toward  the  center  of  the  plasma  column  where  it  is  reflected  to 
form  a  standing  wave. 
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Let  us  turn  now  to  the  profiles  which  model  the  Ormak  plasma. 
This  configuration  is  again  a  Tokamak  plasma  but  as  pointed  out 
above  the  operation  can  be  one  of  two  types,  type  A  and  B 
discharges.   Our  intention  here  is  not  to  discuss  the  nature 
of  these  two  types  of  discharges   but  to  use  only  the  data  result- 
ing from  them.   For  the  Ormak  plasma,  both  A  and  B  type  discharges, 
the  same  results  for  the  cusp  absorption  are  found  as  in  the 
case  of    the  ST  plasma.   The  coil  impedances  are  very  low, 
typically   .01  -  .1  ohms.   The  values   of  impedance  for  the 
Alfven  mode  however  are  very  large   again  and  as  shown  in 
Figure  14,   also  for  a  rather  ideal  case  where  the  coil  now  is 
assumed  to  be  right  at  the  plasma  edge,  can  be  as  much  as  50  ft. 
The  results  shown  in  Figure  14   are  for  a  type  A  discharge  and 
the  wall   has  again  been  taken  to  infinity.   We  notice  some 
difference  from  the  ST  plasma  in  that  the  form  of  the  curves 
for  m  =  1,2,3   show  more  qualitative  differences.   The  calcula- 
tions were  not  carried  far  enough  to  obtain  the  maximum  for  m  =  2 
and  m  =  3.   The  normalized  wave  number,  k,  for  the  coil  was  fixed 
at  .287.   No  interaction  with  stable  kink  modes  was  found  for 
this  case.   If  we  now  consider  the  type  B  profiles  with  all  other 
parameters  held  fixed,  Figure  15,  we  obtain  once  again  strong 
evidence  that  the  spatial  variation  of  the   plasma  profiles 
plays  an  important  role  in  the  actual  energy  absorption.  The 
effect  of  changing  the  coil  wave  length  may  be  seen  in  Figure  16 
where  the  same  parameters  and  profiles  of  Figure  15  were  assumed 
but  with  a  coil  of  exactly  one-half  the  wavelength.   Here  it  can 
be  seen  that  the  wave  length  also  plays  an  important  role  in 
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the  energy  absorption.  Again  no  effect  was  made  to  obtain  an 
optimum  in  the  absorption  rate  with  respect  to  the  coil  wave 
length. 

At  this  point  it  is  appropriate  to  discuss  the  effect  of  an 
externally  conducting  wall.   Here,  we  will  consider  only  a 
perfectly  conducting  wall.   The  position  of  the  wall  may  be  an 
important  consideration   for  the  present  day  and  next  genera- 
tion Tokamak  experiments  for  the  following  reason.    It  is 
desirable  to  create  a  plasma  with  the  largest  possible  cross- 
section  inside  magnetic  field  coils;  the  plasma  column  must  of 
course  be  surrounded  by  a  vacuum  vessel.   The  space  between  the 
edge  of  the  plasma  column  and  the  vacuum  vessel  may  thus  be  very 
small.   Consequently  any  coil  placed  between  the  plasma  column 
and  the  outer  liner  or  vacuum  vessel  could  suffer  from  the 
effects  of  mirror  currents  induced  in  the  outer  wall.  Of  course 
if  the  outer  wall  were  constructed  in  a  segmented  manner 
preventing  long  current  paths,  this  problem  would  be  minimized. 
We  have  performed  a  series  of  calculations  using  the  Ormak  "A" 
profiles  in  which  it  has  been  assumed  that  a  conducting  wall 
surrounds  the  coil  and  have  varied  the  wall  position  from  a 
distance  of  say  10  cm  to  a  few  cm  away  from  the  coil.  Figure  17a 
shows  the  results  from  these  calculations;  shown  is  the  real  impe- 
dance versus  frequency  for  the  case  where  the  coil  is  assumed 
to  be  placed  at  the  edge  of  the  plasma  and  the  wall  position 
varied  from  r  /r   =  10  to  r  /r   =1.1.   As  expected  the  impedance 
decreases  drastically  as  the  wall  is  brought  near  the  coil. 
Noticeable  also  is  the  effect  of  shifting  the  position  of  maximum 
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absorption  to  higher  frequencies  the  closer  the  wall  is  to  the 
coil.   Shown  in  Figure  17b   are  the  values  of  the  imaginary  part 
of  the  coil  impedance  from  which  along  with  the  real  part  of 
impedance,  the  effective  "Q"  of  the  plasma-coil-wall  system 
may  be  calculated.   It  can  be  seen  that  by  dividing  the  imaginary 
by  the  real  impedance,  "Q"  values  of   from  3-6  in  the  low 
frequency  range  are  obtained.   Such  values  allow  the  classifi- 
cation of  such  a  system  as  a  "low-Q"  system  suggesting  efficient 
use  of  energy  in  the  external  circuit.   If  we  now  move  the  coil 
away  from  the  edge  of  the  plasma  and  simultaneously  vary  the 
wall  position  we  obtain  the  results  shown  in  Figure  18.   Here 
the  position   of  the  coil   r  /r   =1.5,  and  again  the  wall 
position  is  varied  from  r  /r   =  10  to  r/r   =  1.1.   Shown  in 
Figure  18a   is  the  real  impedance  versus  frequency  and  character- 
istic decrease  of  impedance  as   the  wall  is  brought  closer  to 
the  coil.   For  this  case  it  is  seen  that  the  position  in  frequency 
of  maximum  absorption  varies  only  slightly   as  the  wall  is  brought 
closer  to  the  coil.   Figure  18b   shows  the  corresponding  computed 
values  of  the  imaginary  part  of  the  impedance  and  again  low  "Q" 
values  are  obtained.   It  is  thus  recognized  that  it  is  desirable 
to  have  ample  space  between  coil  and  wall  in  order  to  insure 
reasonable  power  absorption.   Looking  towards  the  future  reactor 
geometries  it  can  be  seen  that  the  problem  of  the  external  wall 
can  be   minimized  or  perhaps  even  neglected.   This  is  due  to  the 
possibility  of  the  use  of  a  diverter   which  would  necessitate  a 
larger  external  vacuum  region  in  which  coil  structures  could  be 
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placed.   These  considerations  are  not  the  subject  of  the  present 
paper,  and  we  leave  ttiem  to  later  investigations. 

Finally  let  us  consider  the  results  for  the  high-S  configura- 
tion described  earlier.   We  consider  the  6  value  on  plasma  axis 
to  be  .5   and  consider  an  m  =  1   coil  structure.  This  choice 
would  correspond  to  an  I   =  1   system  in  the  actual  helical 
Scyllac  configuration.   The  wave  length  for  the  coil  has  been 
taken  to  be  equal  to  the  helical  wave  length  in  the  present 
Scyllac  device,  and  the  coil  position  r _/r  =1.5.   We  neglect 
here  the  effect  of  the  external  wall.   The  results  of  the  calcula- 
tion of  the  real  impedance  versus  frequency  are  shown  in  Figure  19 
Here,  the  cusp  absorption  rate  shows  a  marked  increase  over  and 
above  that  for  the  previous  Tokamak  configurations  with  maximum 
values  of  the  order  of  one  ohm  or  so.   The  Alfven  absorption  rate, 
however,  again  dominates  as  seen  in  the  figure  and  leads  to  very 
impressive  predictions  concerning  the  coupling  of  energy  to  the 
plasma.   Here  also,  a  calculation  of  the  effective  "Q"  value 
for  the  system  defines  the  process  to  be  a  "low-Q"  system  with 
values  typically  between  3-6  at  maximum  absorption. 

We  have  thus  demonstrated  numerically  that  effective  coupling 
of  energy  from  an  external  source  can  be  accomplished  through  the 
MHD  resonances.   Coupling  is  possible  in  both  the  Alfven  mode 
dominating  for  typical  Tokamak  and  Scyllac  plasma  profiles. 
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10.   Experimental  Verification  of  Alfven  Wave  Heating 

In  this  section,  we  consider  recent  and  past  experimental 
investigations  aimed  at  confirming  various  aspects  of  the  theory 
of  Alfven  wave  damping,  absorption  and  heating.   It  will  be  shown 
that  there  is  ample  evidence  to  confirm  the  theory  set  forth  in 
this  report;  there  exists  close  quantitative  comparison  between 
theoretical  prediction  and  experimental  measurements.   In  order, 
we  will  outline  results  from  the  ISAR  I  helical  theta-pinch 
experiment  (Garching) ,  the  Proto  Cleo  Stellarator  experiment 
(Wisconsin) ,  the  Heliotron  D.  Experiment  (Kyoto) ,  and  a  very 
recent  linear  theta-pinch  experiment  (Lausanne) . 

The  first  full  experimental  verification  which  yielded  a 
confirmation  of  the  fundamental  principles  upon  which  the  theory 
of  resonant  Alfven  wave  heating  is  based  was  performed  in  the 
ISAR  I  helical  theta-pinch  at  Garching.   This  experiment  consisted 
of  a  5.4  meter  long  theta-pinch  coil  whose  inner  surface  was 
machined  to  produce  an  I   =  1  helical  magnetic  field  component. 
The  implosion  phase  of  the  plasma  column  formation  followed  a 
more  or  less  normal  theta-pinch  like  form  with  the  helical 
perturbation  being  felt  at  a  somewhat  later  time.   This  delay 
resulted  in  a  standing  i    =  1  helical  oscillation  of  the  plasma 
column.   The  mode  structure  of  the  column  could  be  described  as 
basically  an  m  =  1  long  wavelength  perturbation.   The  wavenumber 
of  the  perturbation  was  equal  to  that  of  the  helical  field  compo- 
nent,       "•   ,  where  \.  is  the  wavelength  of  the  helical  field. 
H  n 

This  helical  perturbation  was  observed  to  oscillate  at  approximately 
the  Alfven  frequency  u>.    computed  as 
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w. 


(2-6) (kVA) 


(160) 


(the  Alfven  speed  V  was  calculated  from  the  density  inside  the 
plasma  column  and  the  main  magnetic  field  outside  the  plasma 
column).   Further,  the  helical  oscillation  was  observed  to  damp 
away  after  two  to  three  oscillation  peridos.   This  damping  was 
several  orders  of  magnitude  larger  than  that  which  could  be 
expected  to  occur  from  purely  classical  effects  such  as  resist- 
ivity and  viscosity.   The  experimentally  deduced  damping  decrements 
were  compared  with  theoretical  calculations  (Tataronis  and  Grossmann, 
Grossmann  and  Tataronis,  1973)  for  a  wide  range  of  plasma  para- 
meters which  covered  both  collision  dominated  and  collisionless 
regimes.   We  show  below  a  sample  of  the  results  from  the  Garching 
investigation  (Grossmann,  Kaufmann,  Neuhauser,  1973). 
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In  the  above  table,  the  values  of  ion  temperature  T.,  plasma 

density  n  are  given  at  the  time  of  maximum  magnetic  field.   The 

experimental  values  of  y  divided  by  w.  were  approximately  constant 

over  the  entire  regime  from  the  collisional  phase  (cj  .t..  =  .25) 

3  ci  n 

to  the  collisionless  regime  (cj  .x.   =  40)  where  u>   is  the  ion 

Cl  11  ci 

-66- 


cyclotron  frequency  and  T..  is  the  ion-ion  collision  time.   It 
is  noticed  that  the  values  of  (V^A)rrJ.cSICAL  are  many  orders  of 
magnitude  less  than  the  experimental  values.   The  theoretical 
predictions  for  the  damping  assuming  a  simple  form  for  the  in- 
homogeneous  plasma  profiles  can  be  seen  to  agree  well  with  the 
experimental  ones.   These  results  are  felt  to  give  conclusive 
proof  that  the  experimentally  observed  damped  helical  Alfv6n 
oscillations  are  surface  waves  which  decay  due  to  phase  mixing 
(Tataronis  and  Grossmann,  Grossmann  and  Tataronis,  1973). 

The  next  experiment  we  consider  was  performed  in  the  proto- 
cleo  stellarator  at  Wisconsin  (Golovato,  Shohet,  and  Tataronis, 
1976) .   Proto-cleo  is  a  i   =  3  stellarator  with  seven  field  periods 
in  the  toroidal  direction.   The  aspect  ratio  of  the  plasma  is 
between  8  and  9,  a  reasonable  number  for  a  cylindrical  approxi- 
mation to  the  torus.   Alfven  wave  heating  was  produced  with  a 
separate  electrostatically  shielded  helical  winding  whose  pitch 
was  designed  to  approximate  a  q  =  3  rational  magnetic  surface. 
According  to  measured  magnetic  field  and  density  profiles  for  the 
proto-cleo  plasma,  for  a  constant  frequency  of  the  externally 
applied  rf  source,  two  Alfv6n  resonance  regions  were  predicted  at 
minor  plasma  radii  of  approximately  1.5  and  3.5  cm.   During  the 
course  of  the  experiment,  the  resonant  regions  vary  in  position 
according  to  changes  in  plasma  density  profiles.   It  was  experi- 
mentally shown  that  the  incident  Alfven  wave  amplitude  did  not 
decrease  in  the  plasma  region;  this  implies  that  the  wave  penetrated 
into  the  plasma.   Shown  below  are  the  results  of  a  local  measure- 
ment of  the  electron  temperature  with  rf  compared  to  that  without 
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rf  as  a  function  of  the  plasma  minor  radius.   It  is  noticed  that 
a  peaking  of  the  ratio  of  T  with  rf  to  T  without  rf  occurred 
near  the  region  of  the  expected  Alfven  wave  resonance. 
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The  variation  of  energy  absorption  and  heating  with  external 
frequency  was  consistent  with  the  calculations  of  heating 
efficiency  for  proto-cleo  according  to  theoretical  predictions 
(Tataronis  and  Grossmann,  1976)  .   The  temperature  of  both  elec- 
trons and  ions  increased  by  a  factor  of  two  during  heating  and 
this  heating  was  global  in  nature  since  the  coil  structure 
surrounded  the  entire  plasma. 

The  next  experiment  we  consider  was  performed  on  the  heliotron 
D  machine  at  Kyoto  (Uo  et  al ,  1976).   This  device  is  basically  an 
1=2   plasma  with  a  large  rotational  transform  and  strong  magnetic 
shear.   The  aspect  ration  for  this  plasma  is,  similar  to  proto- 
cleo,  about  8.   A  local  I   =   2-like  coil  structure  was  used  to 
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couple  external  rf  power  into  the  Alfven  resonances  of  the  plasma 
column.   Similar  results  to  those  of  proto-cleo  were  observed  in 
the  heliotron  D  device.   That  is,  in  the  regions  where  the 
resonances  were  predicted  to  occur  both  ion  and  electron  tempera- 
tures, were  significantly  increased  (ion  and  electron  temperatures 
were  approximately  doubled) .   It  was  further  shown  that  the  coil 
loading  resistance  was  high  implying  that  the  coupling  corresponds 
to  low-Q  and  broad  resonance.   This  is  consistent  with  theoretical 
predictions  (Tataronis  and  Grossmann,  1976).   The  heating  efficiency 
was  estimated  to  be  as  large  as  60%. 

Finally  we  consider  results  from  an  Alfven  wave  heating  experi- 
ment on  a  high-6  linear  theta  pinch  plasma  column  (Pochelin, 
Keller,  1977).   In  this  experiment  a  rather  dramatic  confirmation 
of  Alfven  wave  heating  was  observed.   The  basic  experiment 
consists  of  a  1.4  meter  linear  theta  pinch  plasma  with  mean 
density  ^  1.3  *  10    cm   ,  total  temperature  ^   50  eV,  and  mean 
plasma  beta  ^  0.25.   A  helical  excitation  winding  giving  an  m  =  1 
perturbation  with  approximately  two  full  wavelengths  was  applied  to 
the  plasma  column.   The  application  of  the  helical  perturbation 
field  took  place  approximately  1.5us  after  the  main  magnetic 
field  had  been  turned  on;  this  allowed  enough  time  for  the  natural 
or  magnetoacoustic  oscillations  to  have  damped  away.   It  was 
noticed  in  the  experiment  that  when  the  exciting  coil  was  driven 
continuously,  that  the  plasma  helical  motion  was  less  damped  at 
frequencies  below  the  Alfven  wave  resonance  peak  (corresponding  to 
the  surface  wave  resonance);  above  this  resonance  point,  the 
helical  plasma  motion  was  damped  rapidly.   When  th«  coil  was 
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operated  with  only  two  excitation  periods  and  at  resonance 
the  kink-like  helical  motion  is  more  strongly  absorbed  than  at 
lower  frequency.   After  the  excitation  was  shut  off,  the  unforced 
plasma  motion  showed  nearly  critical  damping  in  complete  agree- 
ment with  theoretical  predictions  (Tataronis  and  Grossmann, 
Grossmann  and  Tataronis,  1973).   From  a  knowledge  of  the  time 
evolution  of  the  plasma  radius  and  diamagnetic  probe  measurements, 
the  mean  value  of  the  total  plasma  temperature  was  deduced.   The 
plasma  temperature  was  raised  on  the  order  of  15  to  20%  during 
the  short  heating  pulse  and  it  was  further  ascertained  that 
approximately  50%  of  the  excitation  power  was  thermalized  in 
the  plasma.   This  number  also  corresponds  to  a  coupling  efficiency 
of  approximately  50%  and  consequently  shows  the  method  to  be  very 
attractive  for  heating.   The  measured  power  absorbed  by  the  plasma 
as  a  function  of  the  coil  frequency  follows  a  curve  which  is 
identical  to  that  shown  in  Fig.  19. 
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11.   Conclusions 

A  complete,  self-consistent  and  rigorous  treatment  of  rf 
energy  absorption  through  the  use  of  the  linearized  MHD  equations 
of  motion  in  inhomogeneous  plasmas  by  means  of  Alfv^n  waves  has 
been  presented.   A  model  system  is  considered  consisting  of  a 
straight  cylindrically  symmetric  plasma  column  surrounded  by 
a  vacuum  region  in  which  an  idealized  coil  is   placed.  The 
plasma  profiles  are  inhomogeneous  in  the  direction  perpendicular 
to  the  cylinder  axis  and  the  entire  plasma-vacuum-coil  system 
is  surrounded  by  a  conducting  wall.   It  has  been  shown  that  rf 
energy  can  be  irreversibly  absorbed  by  the  plasma  when  the 
external  coil  is  oscillated  at   frequencies  which  lie  in  the 
Alfv^n  continuous  frequency  spectrum  given  by  the   nonuniform 
plasma  profiles.   The  resonant  absorption  mechanism  is  discussed 
in  complete  detail  and  the  following  general  results  have  been 
found:   (i)  strong  absorption  is  predicted  at  reasonably  low 
frequencies  independent  of  the  plasma  total  B  ;  (ii)  variations 
in  plasma  profiles   are  refelcted  in  significant  differences 
in  total  absorption  and  in  the  spatial  position  of  maximum 
absorption;  (iii)  the  position  of  the  coil  and  external  wall 
relative  to  the  plasma  column  plays  an  important  role  in  the  rf 
absorption.   Calculations  with  realistic  Tokamak  and  Scyllac 
plasma  profiles  have  been  carried  out  to  illustrate  the  above- 
values  of  the  effective  coil  impedance  as  a  function  of  coil 
and  plasma  parameters  have  been  computed.   From  the  numerical 
results  it  can  be  shown  that  the  Alfv6n  wave  heating  scheme 
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is  a  low  "Q"  technique  (Q  =  ratio  of  circulating  to  absorbed 
power)  and  should  because  of  the  inherent  low  frequencies  and 
long  wave  lengths  be  considered  for  a  supplementary  heating 
technique  on  future  large  plasma  containment  devices. 

It  is  worth  noting  that  recent  calculations  on 
resistive  (Kappraff,  Tataronis,  Grossmann,  1975;  Kappraff,and 
Tataronis,  1977)  and  kinetic  (Hasegawa  and  Chen,  1975;  Gould, 
1975;  Thompson,  1975)  effects  have  indicated  that  the  energy 
transferred  to  the  resonant  Alfven  region  in  the  plasma  is 
converted  rapidly  to  electron  thermal  energy. 

Finally,  in  support  of  the  above  remarks,  we  mention  that 
experimental  evidence,  in  particular   from  the  proto-cleo  and 
Lausanne  8-pinch  experiments,  has  been  obtained  to  substantiate 
the  claim  that  rf  energy  can  be  coupled  efficiently  to  an  in- 
homogeneous  plasma  column  through  the  continuous  spectrum  of 
stable  Alfven  plasma  oscillations  and  that  the  absorbed  energy 
is  rapidly  converted  to  plasma  thermal  energy. 
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Appendix  A 

The  solution  to  the  MHD  equations  in  the  plasma  regions  is 
matched  to  the  vacuum  electric  and  magnetic  fields  through  eqs. 
(22) -(23).   In  this  appendix  we  perform  this  matching  and  derive 
an  expression  for  the  vacuum  potential  produced  by  the  externally 
applied  vacuum  forces. 

We  express  the  perturbed  vacuum  magnetic  field  5B   as  V<J),  where 
<J>  is  the  magnetic  potential  which  satisfies  Laplace's  equation  in 
the  vacuum  region.   After  Fourier  and  Laplace  transforming 
according  to  eq.  (39),  we  write  the  solution  for  $  in  two  parts 
as  follows, 

$  =  c.  I  (|k|r)  +  c_  K  (|k|r)   ,   r   <  r  •:  r    ,       (Al) 
1   m  2   m  p--c 

$  =  d.  I  (|k|r)  +  d,  K  (|k|r)   ,   r  -:  r  <  r   ,      (A2) 
lm'1       2m'1         c-    -w 

where  I   and  K   are  the  modified  Bessel  functions  of  first  and 
m      m 

second  kind  respectively,  with  argument  |k|r  and  m  interpreted 
in  the  sense  |m|.   The  constants  c  ,  c_,  d1  and  d_  are  determined 
by  the  boundary  conditions  atr=r,r=r,r=r,  namely  eqs. 
(22)-(26).   Application  of  eqs.  (24)-(26)  yields  the  following 
set  of  equations  for  c.,  d.,  and  d_  in  terms  of  c2  and  J  .: 

dl  +  Bd2  =  °  (A3) 

C,     I' (|k|r    )-d,[K' (|k|r    )-   a    I'(|k|r    )]    -   -   C-    K' ( I k | r    )  (A4) 

lmczm  c  mc  m  c 

d.[K  (|k|r  )-  a  I  (|k|r  )  ]  -  c,  I  (|k|r  )  -  c.  K  (|k|r  )  +  — rc-^ 
2   m  '  '  c       m  '  'c       lm''c      2m''c       lk 

(A5) 
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where 


a  = 


K' ( |k  r  ) 
m  '  '  w' 

V  (  |k|r  ) 


(A6; 


and  the  prime,  (  ) ' ,  on  the  Bessel  functions  refers  to  deriva- 
tives with  respect  to  the  argument  |k|r.   Note  that  a  •*■   0  as 

r 


w 


Let  us  define  the  quantity  q  as  follows, 


q  -  ar 


dr 


(A7! 


r=r 


Using  eq.  (A7) ,  the  constant  c~  can  then  be  expressed  in  terms 
of  c,  and  q,  and,  after  solving  eqs.  (A4)  and  (A5)  for  c,  and 
d_  and  rearranging,  the  final  expression  for  <J>  can  be  written  as 


%   =   7— | [I  (Iklr)  -  XiK  (Iklr)]   ,  r   :  r  ■:  r 

1-ax,    m  '  '      Al  m  '  '         p  -    -   c 


where, 


Im(lklrp)-lklq  i;dklrp) 
Km(|k|rp)-|k|q  K;(|k|rp) 


m 


and, 


f  =  -  sgn(k)  i  u0  JsQ  rc  K;<|k|rc) 


1-a 


I' 

m 

(Ik 

^1 

K' 

m 

(|k 

rc}J 

Equations  (A8)  and  (A10)  yield  eqs.  (62)  to  (64) 


(A8) 


(A9) 


(A10) 
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Figure  Captions 

Fig.  1.   Schematic  of  the  plasma-vacuum-coil-wall  configuration. 

Fig.  2.   Friedrichs'  diagram  for  MHD  waves  in  an  infinite  uni- 
form plasma. 

Fig.  3.   Wave  front  diagram  for  the  three  MHD  waves  whose  phase 
velocities  are  shown  in  Fig.  2. 

Fig.  4.   Dispersion  diagram  for  the  three  MHD  waves  of  Figs.  2 

and  3  illustrating  resonance  and  cut-off  features  in  an 
infinite  uniform  plasma. 

Fig.  5.   Schematic  of  the  function  r n(w)  shown  as  a  function  of 
frequency  cj  for  both  the  Alfv^n  and  cusp  branches. 

Fig.  6.   Typical  branch  cuts  and  branch  points  of  the  function 

ty    (r,co)  in  the  complex  oo  plane  (b)  for  the  resonant 

2 
frequency  profile  oj.(r')  (jeA  for  Alfve"n  mode  and  j  =C 

for  cusp  mode)  shown  in  (a) . 

Fig.  7.   Deformed  Laplace  contour  used  for  inverting  the  Laplace 

.  rt      ft 

transform  of  r£(r,w).  Exposed  roots  w  ,  u   in  the  complex 

plane  correspond  to  damped  surface  waves  while  those  on 

the  w   axis  u_,  u  ,  uj   correspond  to  stable  magneto- 
r       m   m   m 

acoustic  waves. 

Fig.  8.   Radial  dependence  of  pressure,  P,  magnetic  fields,  B  /u,  B 

6      z 

and  density  p  in  normalized  form  for  the  ST  Tokamak  plasma. 
Fig.  9.   Radial  dependence  of  Alfven  and  cusp,  ui   and  u>   respec- 

tively,  frequency  continuum  for  the  plasma  profiles  of 

Fig.  8.   The  cusp  frequencies  are  multiplied  by  a  factor 

of  10.   In  the  figure,  m  =  1. 
Fig.  10.  Radial  dependence  of  pressure,  magnetic  field  and  density 

in  normalized  form  for  a  Scyllac  plasma,  3  -  0.5  on  axis. 

-75- 


Fig.  11.  Radial  dependence  of  Alfv£n  and  cusp  continuum  frequencies, 

w-  and  o)r  ,  for  Scyllac  plasma.   The  8  on  axis  equals 

0.5  and  m  =  1  has  been  assumed. 
Fig.  12   Radial  dependence  of  Alfven  and  cusp  continuum  frequencies, 

03-  and  oj   ,  for  reversed  bias  9-pinch.       3  =  0.5 

and      m  =  1  has  been  assumed. 
Fig.  13   Real  part  of  coil  impedance,  Z   ,  versus  frequency, 

f,  for  ST  plasma.   Results  for  m  =  1  and  2  coils  shown 

for  wave  length  equal  to  60  cm. 
Fig.  14   Real  part  of  coil  impedance,  Z   ,  as  a  function  of 

R 

frequency,  f,  for  Ormak  A  discharge.  k  =  .287  and 

m  =  1 ,  2  and  3 . 

Fig.  15   Real  part  of  coil  impedance,  Z   ,  as  a  function  of 

frequency,  f,  for  Ormak  B  discharge.  k  =  .287  and 

m  =  1 ,  2 ,  and  3 . 

Fig.  16   Real  part  of  coil  impedance,  Z   ,  as  a  function  of 

frequency,  f,  for  Ormak  B  discharge,  k  =  .594  and  m=l,2,3. 


Fig.  17a  Real  part  of  coil  impedance,  Z   ,  as  a  function  of  frequency, 

R 

f,  for  Ormak  A  profiles  and  variable  wall  position. 

m  -  1.  k  ■  .287  and  r  /r  varied  from  10  to  1.1.  The 
'  w'  c 

value  of  r  /r   =  1.0  . 
c/  p 

Fig.  17b  Imaginary  part  of  coil  impedance,  ZT  ,  as  a  function  of 

frequency,  f,  for  Ormak  A  profiles  and  variable  wall  position, 
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Fig.  18a  Real  part  of  coil  impedance,  Z_  ,  as  a  function  of 

frequency,  f,  for  Ormak  A  profiles  and  variable 

wall  position,   m  =  1,  k  =  .287  and   rw/r   varied 

from  10   to  1.1.   The  value  of   r  /r   =1.5. 

c   p 

Fig.  18b  Imaginary  part  of  coil  impedance,  Z.  ,  as  a  function 
of  frequency,  f,  for  Ormak  A  profiles  and  variable 
wall  position. 

Fig.  19   Real  part  of  impedance,  Z   ,  versus  frequency,  f,  for 

Scyllac  profiles,   m  =  1,  wave  length  =  33  cm, 

B  =  0.5   and  r  /r   =1.5. 
c     P 
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Figure    11 
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Figure    13 


30  t 


20 


en 

X 

o 

ac 
M 


10" 


4.7 


f  x  I0~6(CPS) 

-92- 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Administration, 
nor  any  person  acting  on  behalf  of  the 
Administration : 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  in  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
in  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
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report. 
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any  information  pursuant  to  his  employment  or 
contract  with  the  Administration,  or  his 
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